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I .  INTRODUCTION 

Characterization  of  the  viscosity  of  hydrocarbon  fuel  mixtures  is 
of  considerable  importance,  because  the  viscosity  directly  controls 
atomization  in  combustors  and  affects  the  performance  of  other  engine 
components  such  as  fuel  pumps.  Typical  aircraft  fuels  are  blends  and 
are  exposed  to  operating  temperatures  which  may  range  from  -55  to  45°C, 
a  temperature  range  which  can  correspond  to  viscosity  changes  of  from 
1  to  3  orders  of  magnitude,  depending  on  the  molecular  weight  of  the  fuel. 
Consequently  an  understanding  of  hydrocarbon  fuel  viscosities  useful  from 
an  engineering  standpoint  demands  an  account  of  the  combined  effects  of 
temperature  and  composition  on  the  viscosity. 

A  truly  fundamental  theory  would  predict  the  viscosity,  along  with 
other  thermodynamic  and  transport  properties,  from  knowledge  of  the  inter- 
molecular  forces  and  radial  distribution  functions  alone.  Such  a  theory 
has  had  appreciable  success  In  application  to  pure,  simple  liquids  such 
as  the  liquefied  rare  gases.  For  solutions,  however,  although  the  general 
theory  has  been  formulated  (1),  it  has  not  yet  been  reduced  successfully 
to  numerical  results.  One  Is  thus  forced  to  approximate  approaches,  of 
which  two  general  types  may  be  distinguished.  The  first  in  which  the 
discrete  molecular  nature  of  the  solvent  is  neglected  has  been  remarkably 


2 


successful  in  explaining  the  viscosity  of  dilute  solutions  of 
high  polymers.  The  second  general  approach  is  to  correlate 
the  viscosity  of  the  mixture  with  the  properties  of  the  pure 
components  and  thermodynamic  parameters  characteristic 
of  the  interaction  between  components. 

Prediction  of  the  viscosity  of  liquid  mixtures  has 
been  a  goal  of  long  standing.  However,  studies  reported 
in  the  literature  on  viscosities  of  organic  and  hydrocarbon 
mixtures  have  never  looked  at  the  temperature  dependence  of 
additivity  rules  over  the  -50  to  50 °C  range  important  for 
fules.  Consequently  we  have  carried  out  a  study  of  the 
effects  of  temperature  and  composition  on  the  shear  viscosities 
of  blends  of  hydrocarbons  of  the  type  found  in  jet  fuel 
kerosenes.  The  combined  temperature-composition  dependence 
of  viscosities  of  five  pure  Cm  hydrocarbons  and  their  different 
mixtures  and  of  binary  Cio  -  Cm  hydrocarbon  system  have 
been  studied,  and  the  results  evaluated. * 

II.  REVIEW  OF  LITERATURE 


The  two  aspects  of  liquid  viscosity  of  interest  here 
have  usually  been  considered  separately  in  the  literature, 
namely,  viscosity  as  a  function  of  temperature  for  pure 
liquids,  and  the  isothermal  dependence  of  viscosity  on 
composition  for  mixtures. 

A.  Simple  Theories  for  Viscosity  as  a  Function  of  Temperature 
for  Pure  Liquids 

*  The  experimental  viscosity  measurement  on  the  binary 

Cio-  Cm  mixture  (i.e.,  mixtures  of  HNN  and  XTHCDPD)  were 
supported  by  Contract  No.  F33615-78-C-2000  from  the  Air 
Force  Systems  Command.  They  are  reported  here  for  the  sake 
of  completeness. 
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A  correct  molecular  theory  of  liquid  viscosity  would  be  based  on 
a  detailed  consideration  of  the  characteristic  molecular  dynamics  of  the 
liquid  state.  Present  theories  of  liquid  viscosity  leave  a  lot  to  be 
desired  in  this  respect  as  has  been  pointed  out  in  the  review  by  Brush  (2). 

The  fundamental  question  that  must  be  addressed  is  how  momentum 
is  transferred  from  one  molecule  to  another.  In  gases  momentum  is  trans¬ 
ferred  primarily  by  molecular  collisions  and  to  a  lesser  extent  by  inter¬ 
atomic  forces.  In  liquids  it  is  an  open  question  as  to  whether  or  not  it 
is  the  attractive  or  the  repulsive  forces  that  are  most  important  in 
controlling  the  viscosity.  What  is  agreed  upon  and  what  is  firmly 
established  in  the  theories  of  Brownian  motion  of  colloidal  particles  is 
that  liquid  viscosity  is  intimately  associated  with  the  diffusion  of  a 
colloidal  particle.  The  diffusion  coefficient  D  and  the  viscosity  n  are 
inversely  proportional  to  one  another.  The  Stokes  -  Einstein  equation 
obtained  from  hydrodynamics  and  the  theory  of  Brownian  motion  adequately 
correlates  experimental  liquid  shear  viscosities  n  and  self-diffusion 
coefficient  D  for  many  liquids  (3a,  3b) 

Dn/kT  =  (1/2tt)  (Na/V) 1  /3  (1) 

where  Na  is  Avogadro's  number,  V  the  molar  volume,  T  the  temperature, 
and  k  the  Boltzmann  constant.  Therefore,  the  problem  of  determining 
momentum  transfer  in  a  liquid  can  be  converted  into  the  conceptually 
simpler  problem  of  determining  thefactors  that  govern  molecular  diffusion. 
Hence  many  viscosity  theories  for  liquids  are  essentially  theories  of 
diffusion.  One  obvious  goal  of  such  theories  is  to  obtain  a  relationship 
between  equatlon-of-state  properties  (pressure-volume-temperature)  and 
the  viscosity  of  a  liquid.  This  makes  it  possible  to  infer  viscosity 


behavior  from  PVT  behavior.  In  this  respect  one  of  the  early  equations, 
and  certainly  one  of  the  most  popular,  relating  the  shear  viscosity  of  a 
liquid  to  temperature  (4,  5,  6)  is  the  so-called  Andrade  equation: 

n  =  A  exp(B/T)  (la) 

where  A  and  B  are  constants.  A  relationship  of  this  kind  was  suggested 
empirically  by  Arrhenius  (5)  and  by  Guzman  (6)  and  derived  theoretically 
by  Andrade  (4). 

Eyring  (7)  gave  the  equation  added  theoretical  significance  in 
the  following  form: 

n  =  h(Xi/X2x,X?)exp(AG*/RT)  (2) 

where  \\ ,  x2 ,  x3  and  x  are  characteristic  distances  shown  in  Fig.  1,  R  is 
the  ideal  gas  constant,  &G*  the  molar  activation  free  energy,  and  h  is 
Planck's  constant. 

If  x  is  assumed  equal  to  Xj,  and  XiX2x3  is  identified  with 
molecular  volume,  Eq.  (2)  may  be  written  as: 

n  =  (hNa/V)exp(6G*/RT)  (3) 

For  kinematic  viscosity  v  Eq.  (3)  can  be  written  as: 

v  =  n/p  =  (hNa/M)exp(AG*/RT)  (4) 

where  M  is  the  molecular  weight. 

Eyring's  picture  of  shear  between  two  layers  of  liquid  involves  the 
successive  passage  of  individual  molecules  from  one  equilibrium  position 
to  another  as  indicated  in  Fig.  1.  Such  a  passage  requires  that  either 
a  hole  or  a  site  be  available.  The  production  of  such  a  site  requires 
the  expenditure  of  energy  to  push  back  other  molecules.  The  movement  of 
the  molecule  may  be  regarded  as  the  passage  of  the  system  over  a  free 
energy  barrier  of  height  aG*. 
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Further,  since  aG*,  the  Gibbs  free  energy  of  activation,  may  be 
replaced  by  AH*  -  TaS*,  it  is  seen  that 

n  =  [(hNa/V)exp(-AS*/R)]exp(AH*/RT)  (4a) 

where  AH*  and  AS*  are  enthalpy  and  entropy  of  activation.  Since  the 
molar  volume  of  a  liquid  V  does  not  vary  greatly  with  the  temperature,  if 
aH*  and  AS*  are  taken  as  constant,  Eq.  (4a)  takes  the  form  of  Eq.  (la). 

Eq.  (la)  predicts  a  linear  plot  of  In  n  vs.  1/T  for  liquids. 
Experimentally,  viscosity  data  which  cover  a  sufficient  range  (^1  order 
of  magnitude  in  ■<  do  not  give  such  plots.  Rather,  In  n  versus  1/T 
plots  are  usually  concave  upwards. 

Some  other  investigators  (Lennard  -  Jones  and  Devonshire  (8); 
Glasstone,  Laidler  and  Eyring  (7);  Frenkel  (9);  Fowler  and  Guggenheim 
(10)  )  have  employed  the  concept  of  the  "Free  Volume"  in  statistical 
thermodynamic  theories  of  the  liquid  state.  According  to  Glasstone, 
Laidler  and  Eyring  the  free  volume  may  be  regarded  as  the  volume  in  which 
each  molecule  of  a  liquid  moves  in  an  average  potential  field  due  to  its 
neighbours.  However,  theoretical  estimates  of  free  volume  depend  on 
postulates  regarding  the  compressibilities  of  the  molecules  and  the  nature 
of  their  packing  in  the  liquid  state.  A  definition  of  free  volume  Vf 
often  used  is  that  employed  by  Doolittle  (11): 

Vf  =  V  -  V0  (5) 

where  V  is  the  measured  volume  at  temperature  T  and  V0, which  is  termed 
the  "occupied  volume", is  the  close  packed  volume  of  the  liquid.  Generally 
Vf  is  presumed  to  vanish  on  isobaric  cooling  at  some  finite  temperature 
To. 

The  free  volume  Vf  is  a  time-average  quantity  and  an  equilibrium 
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property  of  the  system.  In  a  liquid-like  system,  however,  the  local  free 
volume  is  continually  being  redistributed  throughout  the  medium,  the 
redistribution  occuring  simultaneously  with  the  random  thermal  motions  of 
the  molecules.  The  basic  idea  underlying  the  free-volume  approach  is 
that  the  molecular  mobility  at  any  temperature  is  dependent  on  the  avail¬ 
able  free  volume  at  that  temperature.  As  temperature  increases,  the  free 
volume  increases  and  molecular  motions  become  more  rapid. 

Doolittle  (11)  found  that  the  temperature  dependence  of  the  shear 
viscosity  of  low  molecular  weight  hydrocarbon  liquids  could  be  represented 
by  an  empirical  equation  of  the  following  form: 

n  =  a  exp  (b/f)  (6) 

where  a  and  b  are  constants  and  f  =  Vf/V  is  the  fractional  free  volume. 

Cohen  and  Turnbull  (12)  gave  this  expression  a  theoretical  basis 
by  showing  that  Eq.  (6)  can  be  derived  by  assuming  a  distribution  of  hole 
sizes,  so  that  the  jump  probability  is  determined  only  by  the  chance  of 
a  molecule  finding  an  adjacent  local  free  volume  of  sufficient  size  to 
jump  into.  Assuming  a  quasicrystalline  liquid  structure,  a  molecule  can 
be  pictured  as  vibrating  about  an  equilibrium  position  until  a  combina¬ 
tion  of  two  events  occurs:  (1)  the  molecule  attains  sufficient  energy 
to  overcome  the  attractive  forces  holding  it  to  its  neighbors,  and  (2) 
an  empty  site  is  available  into  which  the  molecule  can  jump.  Therefore, 
the  probability  of  a  transition  or  jump  from  one  site  to  another  Pj  is 
given  by  the  expression 

Pj  =  PE  x  Py 

where  P^  is  the  probability  of  attaining  sufficient  energy  to  break  bonds 
and  Pv  is  the  probability  that  there  is  sufficient  local  free  volume  for 
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a  jump  to  occur.  Cohen  and  Turnbull  consider  a  liquid  consisting  of 
hard  spheres  with  only  repulsive  forces.  Each  molecule  is  confined  in  a 
cage  made  up  of  its  neighbors.  For  diffusive  motion  to  occur  there  must 
appear  a  fluctuation  in  the  local  density  which  opens  up  a  hole  large 
enough  to  permit  a  significant  displacement  of  the  molecule  within  the 
cage.  They  assume  such  fluctuations  occur  without  energy  change,  i.e., 

P£  =  1 .  They  associate  a  local  free  volume  with  each  molecule  and 
assume  a  distribution  of  local  free  volumes  or  hole  sizes.  In  this  way 
the  molecules  can  be  characterized  by  their  associated  free  volumes. 

They  next  calculate  an  expression  for  the  excess  entropy  due  to  this 
distribution  of  local  free  volume  V  and  then  find  the  distribution  of 
v  by  maximizing  this  excess  entropy  holding  the  number  of  molecules  and 
total  free  volume  constant.  They  find  the  distribution  of  hole  sizes 
P(v)  t0  be 

P(v)  =  (y/vf)exp(-Yv/vf )  (7) 

where  P(vj  is  the  probability  of  finding  the  free  volume  v  nearby.  The 
average  free  volume  per  molecule  is  Vf.  The  constant  y  is  a  numerical 
factor  needed  to  correct  for  the  overlap  of  free  volume  between  molecules 
and  lies  between  1/2  and  1.'  Assuming  that  some  minimum  local  free  volume 
v*  is  necessary  for  a  jump  to  occuc  one  can  calculate  the  probability  of 
finding  v*  and  thus  the  jump  probability.  This  is  of  course  simply  the 
quantity  Pv  defined  above  and  is  given  by 

tn 

Pv  =/v*  P(v)dv  =  exp(-yv*/vf)  (8) 

Since  the  quantity  v*  should  be  close  to  vo ,  the  close  packed  molecular 
volume,  the  hard-sphere  viscosity  equation  takes  the  form 


n  =  A/Pv  =  A  exp(yv.>/vf) 


(9) 


is 


Cukierman,  Lane,  and  Uhlmann  (13)  have  adopted  the  formulation  due 
to  Turnbull  and  Cohen,  and  for  the  case  of  small  free  volume,  they  have 
expressed  the  self-diffusion  coefficient  for  transport  in  a  liquid  by 

D  -  (  iuV*/3)exp((-, V*/V| )  (10) 

I  /  ' 

where  u  is  the  gas  kinetic  velocity  and  may  be  expressed  by  u  =  (3kT/m) 
where  m  is  mass  of  the  molecule;  and  a  is  a  geometric  factor  relating  the 
molecular  displacement  to  the  free  volume. 

In  relating  this  expression  to  the  shear  viscosity,  it  is  customary 
to  assume  that  D  is  inversely  proportional  to  u: 

D  =  b/n  (11) 

For  relatively  simple  organic  liquids,  the  Stokes  -  Einstein  expression 
provides  a  useful  estimate  of  b,  at  least  in  the  range  ot  high  fluidity: 

b  =  kT/3.iao  (12) 

where  a«  is  the  molecular  diameter.  Using  Fqs.  (10)  -  (12),  and  taking 
uV  +  pao  and  y V*  =  <s V,n  (Vm  being  the  molecular  volume).  Cukierman,  Lane 
and  Uhlmann  have  expressed  the  viscosity  as 

n-  (0.2/»au‘ )(ilikT)1/-,exp(.*-Vm/Vf)  (13) 

For  a  fractional  free  volume  of  fa  at  temperature  Td,  and  assuming  that 
the  glassy  expansion  coefficient  provides  a  good  estimate  of  the  van  der 
Waals  expansion  in  the  liquid-liquid  range,  Eg.  (13)  becomes 

n  -  (0.2/na./>)(mkT)i/'expia/(fd  i.\„  (T  -  Td)|.  (14) 

This  expression  is  equivalent  to  a  modified  Vogel  -  Fulcher  expression 

loyM,  r,  --  A  ♦  (  1/2) log  T  •  li/(T  -  T.)  (15) 

where  10^  =  (0.2/i  a-/  ) (mk ) * ^‘  ,  li  <■■/?.  ,  and  To  =  Td  -  (fj/. ■.>.), 
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i.e.  To  is  the  temperature  where  tree  volume  vanishes  on  isobaric 
cool iny . 

Recently,  Laugh! i n  and  Uhlmann  (13)  have  reported  a  study  of 
viscous  flow  in  a  number  of  relatively  simple  organic  liquids.  They 
have  concluded  that  for  none  of  the  liquids  investigated  by  them,  the 
flow  behavior  over  the  full  range  of  measured  viscosity  can  be  adequately 
represented  by  any  of  the  standard  theoretical  models.  However,  they 
note  that  agreement  is  found  between  molecular  dynamic  calculations  and 
a  modified  free  volume  theory.  In  addition,  the  free  volume  theory  gives 
a  good  description  of  viscous  flow  at  high  temperature.  Hence  they 
suggest  that  free  volume  theories  are  most  appropriately  used  to  describe 
flow  in  the  high-temperature  region  rather  than  in  the  region  around  the 
glass  transition.  They  have  also  suggested  that  a  second,  easier  flow 
process  becomes  dominant  in  the  low  temperature  region. 


III.  VI SCOS I TY  OF  MI XTURES 

The  viscosity  of  liquid  mixtures  has  attracted  much  attention  in 
the  literature,  both  from  the  practical  and  theoretical  standpoints. 

But  it  should  be  mentioned  that  none  of  those  papers  have  considered 
mixture  viscosities  over  a  sizable  temperature  range. 


1.  "Ideal"  liquations 
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Here  nmi-x,  n (  and  are  ttie  viscosities  of  the  mixture  and  pure 
components  1  and  2  respectively,  X(  and  X,  are  the  mole  fractions  of 
the  components  1  and  2  in  the  mixture.  It  may  be  pointed  out  that  these 
and  similar  other  simpler  forms  have  not  proved  successful  even  for  the 
prediction  of  the  viscosities  of  mixtures  of  yases  at  ordinary  pressures. 

2.  Equations  Based  on  Transition  State  Theory 

a.  One  parameter  equation  for  binary  mixtures 
Katti  and  Chaudhri  (16)  have  assumed  the  transition  state  theory 
expression  for  viscosity  and  assumed  that  the  activation  free  energies  of 
mixtures  can  be  approximated  by  a  regular  solution  theory. 

The  general  form  of  the  transition  state  expression  for  viscosity 
obtained  by  rearranging  Eq.  (3)  is 

aG*/RT  =  ln( .  V/hNa) (■ (19) 

For  the  mixture,  they  have  assumed: 

AG*mix/RT  -  X ,  ( AG  *  /RT)  +  X,(aG*  /RT)  +  X^.JW/RT)  (20) 

where  W  is  the  excess  heat  of  mixing  parameter  that  appears  in  regular 
solution  theory.  Also  they  assume: 


1 1 

ln  Unix/-' mix)  ~  V1'  (“j/-,)  1  X.  In  (.,/a  J  (21) 

Using  Eq.  (19)  for  the  respective  mixture  ami  pure  component  viscosities 
then  (jives 

log  -imixVmix  -  log  n1Vi  *  X  log  -,  V  +  X^JW/RT)  (22) 
where  tlie  V‘s  ore  the  molar  volumes. 

b.  Two  parameter  equation 

Heric  (17)  has  extended  the  regular  solution  approach  of  Katti 
and  Chaudhri  to  multicomponent  systems  and  obtained  tire  following 
expression  for  kinematic  viscosity  ■.> : 

n  n  n 

log  v  =  Xjloy  vj  +  r.  Xjlog  Mj  -  log  ).  X^Mi  +  „  (23) 

i=i  i=i  i=i 

If  only  binary  interaction  terms  are  assumed  necessary,  then  Katti  and 
Chaudhri' s  theory  gives: 

*i  =  (1/2)"  "  X|X|U1J  (24) 

• i=i i=i  J  J 

where  6^  is  a  deviation  function,  representing  departure  from  a 
noninteracting  system,  ajj  is  an  interaction  parameter  with 
and  u{]  =  .*jj  =  0. 

Two  aspects  of  Eq.  (23)  must  now  be  taken  into  account.  The 
first  one  -  the  use  of  a  single  <jj  term  for  representation  of  each 
binary  interaction  -  has  been  considered  by  lleric,  using  the  data  of 
Katti  and  Chaudhri.  lie  found  that  Eq.  (24)  did  not  give  an  adequate  fit 
to  the  data  and  concluded  that  in  a  multicomponent  system,  assuming 
binary  interaction  only,  a  relationship  containing  more  parameters  than 
Eq.  (24)  was  required: 


i .  .  n  - .  x  i  A  jl  '  i  l  '  ‘  i  j  ‘  Ai  '  x  j )  I 

1  j 

i  .1 

1  he  second  aspect  of  t  * | .  { 2A )  -  the  >nlei|ii(U  >  ot  assuming  binary  inter¬ 
actions  only  -  has  been  considered  again  by  llei  ic  (17)  with  viscosity 
dci  t  <i  ori  the  t.euiwry  system  reported  by  Kalidas  and  faddha  (Id).  He 
tound  that  the  as  sump t  ion  of  only  binary  interactions  inadequate  to 
describe  the  data  exactly  and  for  a  ternary  system  Lq.  (?i>)  must  ho 
expanded  by  inclusion  of  a  term  for  ternary  interaction. 

It  may  be  anticipated  that,  for  systems  higher  in  number  of 
components,  higher  interaction  terms  may  well  be  required.  Thus  in 
quaternary  systems,  for  example,  a  tern  in  X,X .  X  X,t  would  be  required. 

3.  Free  Volume  Approach 


It  was  mentioned  previously  that  theories  relate  the  viscosities 
of  liquids  either  to  the  activation  energy  required  for  the  molecule  to 
overcome  the  attractive  forces  of  its  neighbors  and  flow  to  a  new  position 
(absolute  reaction  rate  theory)  or  to  the  probability  that  an  empty  site 
exists  near  a  molecule  (free  volume  theory).  Macedo  and  litovitz  (19) 
have  made  the  hypothesis  that  the  two  effects  are  combined,  so  that  the 
probability  for  viscous  flow  is  taken  as  the  product  of  probability  for 
acquiring  sufficient  activation  energy  and  ot  the  occurrence  ot  an 
empty  site.  Similar  assumptions  can  be  made  for  solutions. 

Upon  combining  the  absolute  reaction  rate  and  free  volume 
theories,  one  obtains  (19)  for  the  viscosity  of  the  solution 


/V,)] 


n 


-  A  exp[ (  r.‘/l<r)  M,v* 


(26) 
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An  equation  of  the  same  to  tin  hold-,  for  the  pure  components 

ni  =  A  exp[( \Gj*/Rl )  f  (yV*/Vfj)J  ( i=I  .2)  (27) 

where  the  subscript  labels  the  property  of  pure  component  i.  Considering 
the  fundamental  assumptions  of  absolute  rate  theory  and  following 
Rosevear,  et  al.,  (20)  one  now  assumes 

AG*  =  Xj.'.G1*  +  X,;.g/  -  ,.aGr  (28) 

where  aGR,  the  excess  free  energy  of  activation,  is  closely  related  to 
the  excess  free  energy  of  mixing,  and  u  is  a  constant  of  order  unity. 

Substitution  of  Eq.  (28)  in  equation  (26),  taking  logarithms  of 
the  resulting  equation  and  of  Eq.  (27),  and  simple  manipulation  then 
yields 

In  =  XJn  Hj  +  X,ln  r,2  -  (..GR/KT)  +  V*C(l/Vf)  -  (X,/Vfj) 

-  (X2/Vf2)]  (29) 

The  excess  free  energy  of  activation  can,  in  turn,  be  broken  down  into 
enthalpy  and  entropy  contributions 

AGR  =  AHM  -  TaSR  (30) 

* 

R 

aHm  is  the  enthalpy  of  mixing  per  mole  of  solution,  and  AS  is  the 
residual  entropy  per  mole. 

The  reduced  volume  could  be  defined  as: 

V  =  V/V*  (31) 

Substitution  of  Eqs.  (5),  (31)  and  (30)  into  Eq.  (29)  yields 

In  r,  =  X  In  n,  +  X.ln  n2  -  (aHm/RT)  +  (aSr/R)  +  l[l/(V  -  1) 

-  Xt/( Vj  -  1)  -  y(V2  -  1)]1  In  n1d  *  In  n,|  +  In  ns+  In  ny  (32) 
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This  equation  displays  explicitly  the  various  contributions  to  the 
mixture  viscosity:  the  ideal  viscosity,  tin-  enthalpy  and  residual 
entropy  of  mixing,  and  the  difference  in  free  volume  lie tween  solution 
and  pure  components. 

Jambon  and  Dolmas  (21)  have  obtained  \GR  from  the  experimental 
■.lift  and  the  calculated  \Sf^.  lhe  expression  for  ,‘.Sp|R,  the  excess 
entropy  of  mixing,  is 

aSMR  3X,  (Pj  *Vl  */Rl ,  *)ln[(VJ7  5  - 1 ) / ( V '  7  '  -  1 )  J  •  3X  (P  *V  *) 

/(RT.*)ln[(vy ‘  -1 J/CV1 7  ’  -  1)]  (33) 

P*,  V*,  and  T*  are  respectively  the  reduced  pressure,  volume,  temperature 
reduction  parameters.  They  are  obtained  tor  the  pure  component  through 
the  equations  V  =  (4.J/3  +  1)/(..T  i  1)  and  V*  -  V/V  where  •<  is  the 
expansion  coefficient.  T  is  obtained  from  1  -  V ~ 1 ( 1  -  V-1/1)  and  T*  - 
T/T.  P*  is  obtained  from  the  thermal  pressure  coefficient  , :  P*  = 
yV  T. 

According  to  Bloomfield  and  Dewan  (23),  it  is  not  clear  a 
priori  whether  all  of  the  contributions  to  the  mixture  viscosity  in 
Eg.  (32)  are  equally  important  and  should  be  considered  together  in 
computing  r, .  They  have  shown  that  the  experimental  viscosities  are 
best  reproduced  by  either  the  absolute  rate  theory  alone  or  the  free 
volume  theory  alone,  for  which  agreement  within  theory  and  experiment 
is  generally  within  a  few  percent. 


b 


l V .  EXPERIMENTAL  PROCLDURE 

A.  Materials. 

Seven  liquids  were  used  in  our  study:  six  C10  hydrocarbon 
liquids  -  n-Decane  (ND),  n-Butyl cycl ohexane  (NBC),  ci s-Decaliydronaph- 
thalene  (CD),  n-Butylbenzene  (NBB)  ,  2 ,7-Di methyl octane  (DM0),  exo- 
tetrahydrodicyclopentadiene  (XTHDCPD),  and  one  hydrogenated  dinier  of 
norbornadiene  (hexacyclic  endo-endo-dihydrodinorbornadiene ,  HNN)  The 
structures  of  these  compounds  and  their  formulae  are  given  below. 


ND 

C10H22 


o 


NBB 

C10H14 


HNN 

c14h18 


XTHDCPD 

C10H16 


Samples  of  n-Decane  (99+%),  n -Butyl cyclohexane  (99+%),  cis-Decahydronaph- 
thalene  (99%),  and  n-Butylbenzene  (99+%)  were  obtained  from  the  Aldrich 
Chemical  Company,  Milwaukee,  Wisconsin.  2 ,7-Dimethyl  octane  (99%)  was  ob¬ 
tained  from  Chemical  Samples  Company,  Columbus,  Ohio.  XTHDCPD  (99.9%) 
was  obtained  from  Dr.  A.  Schneider  of  Suntech,  Inc.  These  were  used  with 
out  further  purification.  Impure  HNN  was  obtained  from  Dr.  A.  Schneider 
of  Suntech,  Inc.  Pure  HNN  (99.9%)  was  prepared  from  M.  Shahriari  and 
R.  Mossadegh  by  triple  recrystallization  from  acetone  at  dry  ice  tempera¬ 
ture,  followed  by  vacuum  distillation. 


For  the  first  five  Cl0  hydrocarbons  (ND,  NBC,  CD,  NBB  and  DM0) 
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shear  viscosities  n  and  densities  p  were  measured  for  the  five  neat 
liquids,  all  ten  binary  0.5  -  0.5  mixtures,  four  ternary  1/3  -  1/3  -  1/3 
mixtures,  one  quaternary  0.25  -  0.25  -  0.25  -  0.25  mixture,  and  one 
quinternary  0.2  -  0.2  -  0.2  -  0.2  -  0.2  mixture.  In  addition,  n  and  p 
were  determined  for  the  0.25  -  0.75  and  0.75  -  0.25  CD/DM0  and  NBC/DM0 
binary  mixtures.  Finally  n  and  p  were  determined  for  the  pure  HNN 
XTHDCPD  and  four  different  (0.8  -  0.2,  0.6  -  0.4,  0.4  -  0.6,  and  0.2  - 
0.8)  HNN/XTHDCPD  binary  mixtures.  All  of  the  above  compositions  are  in 
mole  fraction.  The  temperature  ranges  of  measurement  were  from  50°C  down 
to  either  about  -55°C  or  the  temperature  above  this  where  the  sample 
crystallized.  The  melting  points  of  pure  NBC  (-75°C),  DM0  (-55°C), 

NBB  (-88°C) ,  and  XTHDCPD  (-91°C)  are  at  the  lower  limit  of  or  below  the 
temperature  range  of  this  studty.  The  melting  points  of  CD  (-43°C)  and 
of  HNN  (8°C)  are  within  this  range,  but  CD  and  HNN  were  stable  under 
supercooling.  CD  was  never  observed  to  crystallize,  while  measurements 
on  pure  HNN  could  be  extended  down  to  -51 °C  without  crystallization. 
HNN-XTHDCPD  mixtures  did  not  crystallize.  ND  (mp  =  -30°C)  and  most  of 
the  mixtures  containing  substantial  amounts  of  it  could  be  supercooled 
very  little,  if  at  all,  and  crystallized  quite  reproducibly  whenever 
the  temperature  dropped  below  a  certain  critical  value.  These  tempera¬ 
tures  were  about  -31°C  for  pure  ND,  -40°C  for  0.5  ND  -  0.5  NBC,  -44°C 
for  0.5  ND  -  0.5  CD,  -45°C  for  0.5  ND  -  0.5  NBB,  -51°C  for  1/3  ND  - 
1/3  NBC  -  1/3  CD  and  -51°C  for  1/3  ND  -  1/3  CD  -  1/3  NBB. 


B.  Experimental  Apparatus  and  Procedure 
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Kinematic  viscosities v  were  measured  with  factory  calibrated 
Cannon  -  Ubbelohde  semi-micro  capillary  viscometers  from  the  Cannon 
Instrument  Co.,  State  College,  Pennsylvania.  Six  different  viscometers 
of  constants  0.004142,  0.01681,  0.0937,  0.2515,  1.099  and  7.98  cSt/s 
were  used  in  this  study.  To  prevent  the  condensation  of  moisture  inside 
the  viscometers  at  low  temperature,  all  openings  of  the  viscometers 
were  protected  with  small  drying  tubes  containing  anhydrous  calcium 
sulfate  (Drierite). 

For  highly  viscous  liquids  such  as  HNN  and  some  of  its  mixtures 
the  viscometer  efflux  times  were  in  the  range  100  to  1700  s.  For  lower 
viscosity  liquids,  the  efflux  times  were  usually  between  60  and  600  s. 

Each  efflux  time  was  measured  at  least  twice  and,  more  commonly,  at  least 
three  times  with  a  repeatibility  of  0.1  to  0.2%.  Prior  to  each  experiment 
the  viscometers  were  rinsed  with  acetone  to  dissolve  any  hydrocarbon 
remaining  from  a  previous  experiment,  followed  by  cleaning  with  detergent, 
rinsing  with  deionized  water  and  drying  at  110°C. 

Densities  p  were  measured  at  the  same  time  as  the  viscosities  in 
a  dilatometer  of  about  H  ml  volume  constructed  by  sealing  a  portion  of 
a  graduated  5  mL  pipette  stem  to  a  10  mL  volumetric  flask.  Liquid  volumes 
could  be  read  to  the  nearest  0.01  mL.  The  dilatometer  was  calibrated 
using  deionized  water  as  a  density  standard. 

The  viscometers  and  dilatometer  were  thermostatted  to  +0.05°C 
or  better  in  a  stirred  and  regulated  water  bath  above  room  temperature 
and  In  a  stirred  and  regulated  refrigerated  4-liter  methanol  bath 
contained  In  an  unsilvered  dewar  flask  below  room  temperature.  Except 
for  a  narrow  window  to  allow  observation  of  the  dilatometer  and  viscometer, 
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the  methanol  bath  was  fully  covered  with  aluminum  foil  for  better 
prevention  of  heat  leakage. 

Bath  temperatures  were  determined  with  an  accuracy  -  0.02°C 
with  a  copper-constantan  thermocouple  which  had  been  calibrated  by  us 
against  a  Leeds  &  Northrup  model  8167-25  Pt  resistance  thermometer, 
serial  no,  1812215.  The  Pt  thermometer  had  been  calibrated  at  National 
Bureau  of  Standards  at  the  triple  point  of  H20  and  the  tin  and  zinc 
melting  points.  Thermocouple  emf's  were  measured  with  a  Leeds  & 

Northrup  type  K-4  potentiometer. 

Mixtures  with  a  specified  composition  (in  mole  fraction)  were 
prepared  by  direct  weighing  from  the  pure  components. 


Experimental  Results  and  Data  Analysis 


Density 

The  experimental  values  of  density  p  for  all  the  liquids 
studied  are  given  in  Table  I.  The  liquid  densities  p  of  all  the  pure 
compounds  and  their  mixtures  were  found  to  vary  linearly  with  temperature. 
Density  results  are  given  in  Table  II  in  the  form  of  parameters  obtained 
from  least  squares  fits  of  the  experimental  data  to  the  equation 

p(g/Cm3 )  =  p(o)  -  bT(°C)  (34) 

where  p(o)  Is  the  density  at  0°C.  The  volume  measurement  is  the  only 
Important  source  of  error  in  the  density.  The  overall  accuracy  of  the 
density  Is  estimated  to  be  0.1%.  Duplicate  density  results  (Runs  I 
and  II)  were  obtained  for  five  of  the  entries  In  Table  II  and  agree 
within  our  estimated  uncertainty  of  0.1%  or  better,  as  may  be  seen  by 


Comparison  of  results  of  duplicate  viscosity  experiments 


Molar  Composition 

T(°C) 

n(CP) 
Run  1 

Run  2 

%  diff 

NO 

50 

0.612 

0.614 

0.3 

0 

1 .291 

1 .292 

0.1 

-30 

2.572 

2.569 

0.1 

NBC 

50 

0.836 

0.839 

0.4 

0 

1  .889 

1  .890 

0.1 

-30 

4.16 

4.16 

0.0 

-55 

11  .35 

11  .37 

0.2 

CD 

50 

1  .848 

1.844 

0.2 

0 

5.53 

5.51 

0.4 

-30 

15.67 

15.59 

0.5 

-55 

57.1 

56.7 

0.7 

NBB 

30 

0.897 

0.905 

0.9 

0 

1  .469 

1  .467 

0.1 

-30 

3.036 

3.010 

0.9 

-55 

7.69 

7.66 

0.4 

0.5N0  -  0.5NBC 

50 

0.703 

0.700 

0.4 

0 

1.523 

1 .524 

0.1 

-30 

3.147 

3.140 

0.2 

-40 

4.32 

4.29 

0.7 
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comparing  the  values  of  p(o)  and  the  density  temperature  coefficients  b. 

Shear  Viscosities 

The  experimental  values  of  shear  viscosity  n  are  listed  in  Table 
I  for  all  the  liquids  studied.  In  Figs.  2,  2a,  2b  are  shown  the  typical 
examples  of  the  temperature  dependence  of  the  shear  viscosity  in  form  of 
Arrhenius  plots  of  log  n  versus  103/T(K).  The  plots  are  curved,  but 
may  be  described  well  by  the  semi -enipi ri cal  Vogel-Tammann-Fulcher  (VTF) 
equation. 

In  n(P)  «  A  +  B/[T(K)  -  To]  (35) 

where  A,  B,  and  To  are  constants.  The  n  versus  T  data  for  all  liquids 
studied  were  fit  to  Eq.  (35)  in  the  following  fashion.  A  value  of  To 
was  chosen,  the  corresponding  parameters  A  and  B  evaluated  by  a  linear 
least  squares  fit  of  In  n  versus  1/(T  -  To)  and  the  standard  deviation, 
Std.  Dev.  In  o,  of  the  experimental  data  from  the  least  squares  line 
calculated.  The  value  of  To  was  changed  and  the  fit  repeated  until  the 
value  of  To  (to  the  nearest  1  K)  correspond  to  a  minimum  in  Std.  Dev. 

In  n  was  located. 

In  Table  III  the  best  fit  parameters  for  Eq.  (35)  are  listed 
for  the  liquids  whose  viscosity  was  studied.  Std.  Dev.  In  n  is 
typically  about  0.003  for  the  5  C1Q  hydrocarbons  and  their  mixtures, 
and  about  0.004  for  HNN,  XTHDCPD  and  their  mixtures.  Since  Ain  n  = 

An/n,  this  means  that  the  scatter  of  the  data  from  the  best  fit  curves 
is  typically  about  0.3%  and  0.4%,  so  that  the  parameters  of  Table  III 
may  be  used  to  calculate  the  viscosities  of  the  respective  liquids 
rather  precisely  over  the  temperature  ranges  indicated. 

The  shear  viscosities  reported  here  are  accurate  to  within  0.5%  . 
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This  accuracy  has  been  estimated  from  the  corresponding  uncertainties 
in  temperature  control  and  in  the  measurements  on  the  times  of  flow,  the 
densities  and  the  viscometer  calibration  constants. 

In  Table  IV  shear  viscosities  and  densities  of  some  of  the 
liquids  studied  in  this  work  are  compared  with  those  of  a  number  of 
different  investigators.  At  a  given  temperature  the  densities  measured 
in  different  laboratories  agree  with  ours  well  within  the  experimental 
error.  The  disagreement  between  our  viscosity  measurements  and  those 
of  other  investigators  is  at  worst  2%.  This  is  acceptable  in  view  of 
the  many  oossible  sources  of  error  in  viscosity  measurements. 


V .  DISCUSSION  OF  TEMPERATURE  /  COMPOSITION  DEPENDENCE  OF  VISCOSITY 


A.  ^Preliminary  Connients 

We  consider  here  only  simple  equations  with  few  adjustable 
parameters  to  describe  composition  and  combined  composition/temperature 
dependence  of  viscosity. 

B.  "Ideal"  Equations 

We  tested  three  simple  additivity  rules  for  expressing  mixture 
viscosities  In  terms  of  the  viscosities  of  the  pure  components. 

1.  Additivity  of  viscosities: 


f 
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n  =  i .  X  -j  n  i  (36) 

2.  Additivity  of  fluidities: 

1/n  =  EX-j/ni  (37) 

3.  Logarithmic  additivity  of  viscosities  or  fluidities: 

lnn=£Xilnni  (38) 

where  Xi  is  the  mole  fraction  and  ni  the  viscosity  of  pure  component  i. 

Eq.  (36)  grossly  overestimated  mixture  viscosities  and  will  not 
be  discussed  further.  Tests  of  Eqs.  (37)  and  (38)  are  illustrated  in 
Table  V,  where  ratios  of  calculated  to  experimental  mixture  viscosities 
are  tabulated  at  two  temperatures  and  in  Figs,  3,  4,  5,  where 
isotherms  of  log  n  vs.  X^q,  Xrbc  and  XHNN  are  plotted  for  the  CD-DMO, 
NBC-DMO  and  HNN-XTHDCPD  binary  systems. 

Eq.  (38)  uniformly  tends  to  overestimate  mixture  viscosities, 
and  the  errors  become  larger  the  greater  the  difference  between  the 
viscosities  of  the  pure  components,  i.e.,  the  lower  the  temperature. 

Eq.  (37),  which  assumes  isothermal  additivity  of  fluidities, 
gives  the  best  agreement  with  the  experimental  viscosities.  Two 
general  trends  are  immediately  apparent  from  the  ratios  of  calculated 
to  experimental  viscosity  in  Table  V. 

1.  For  binary  mixtures  of  the  saturated  hydrocarbons  (ND,  NBC, 

CD,  DM0)  and  also  for  the  binaries  of  HNN-XTHDCPD,  Eq.  (37)  uniformly 
underestimates  viscosities,  and  the  errors  become  larger  the  larger 
the  differences  In  the  pure  component  viscosities,  i.e.,  the  lower  the 
temperature. 

2.  For  binary  mixtures  of  the  aromatic  hydrocarbon  NBB  with 

the  saturated  hydrocarbons  Eq.  (37)  overestimates  the  mixture  viscosities. 
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These  two  effects,  i.e.,  viscosity  underestimates  which  become  larger 
the  greater  the  difference  in  component  viscosities  and  viscosity 
overestimates  when  a  saturated  and  aromatic  hydrocarbon  are  mixed, 
compensate  one  another  in  mixtures  in  which  both  effects  are  important. 

For  example,  in  the  quaternary  mixture  Eq.  (37)  correctly  predicts  the 
viscosity  within  21  or  better.  This  is  apparently  due  to  the  compensa¬ 
ting  effects  of  the  high  viscosity  component  CD  and  the  aromatic  NBB. 

C .  Katti  and  Chaudhri  Approach  for  Binary  Mixtures 

The  Katti  and  Chaudhri  model  (16)  should  be  capable  of  accounting 
for  the  combined  temperature/composition  dependence  of  binary  mixture 
viscosities  using  a  single  adjustable  parameter,  W.  We  can  rewrite 
their  expression  (Eq.  (22))  in  the  form: 

W/R  =  [T/yi  -  X2)-J[ln  Y  -  Xj  In  Yj  -  X2ln  Y2  ] 

where  Y  stands  for  either  n  or  nV.  According  to  Katti  and  Chaudhri's 
theory,  W/R  should  be  a  constant  over  the  full  range  of  temperature  for 
any  binary  mixture. 

We  have  calculated  values  of  W/R  from  the  experimental  viscosity 
data  for  three  binaries.  The  results  are  shown  in  the  Tables  VI,  VII, 
and  VIII.  To  estimate  the  uncertainty  in  W/R,  we  have  used  the  expression, 

[T/X2(l  -  x2)][o2  +  XlV  +  x2°2231/2 

where  a  is  the  uncertainty  in  In  Y  and  similarly  for  o,  and  a? .  o,  a, 
and  o2  were  set  equal  to  the  standard  deviations  from  the  best  fit  expres¬ 
sions  to  the  VTF  equation  (Table  III). 

For  the  two  binaries  (CD-DMO  and  HNN-XTHDCPD)  in  which  there  are 
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large  differences  in  the  viscosities  of  the  two  pure  components  and  in 
which  there  are  large  deviations  from  additivity  in  In  n ,  W/R  is  not 
constant  as  predicted  by  this  model,  but  shows  definite  trends  within 
experimental  error.  In  particular,  at  constant  temperature  W/R 
increases  with  increasing  mole  fraction  of  the  more  viscous  component. 
At  constant  composition,  W/R  increases  with  decreasing  temperature. 

In  n  in  the  NBC-DMO  system  is  nearly  additive,  so  that  the  W/R 
values  are  small.  Hence  in  this  system  no  trend  in  W/R  can  be  detected 
within  experimental  error. 

D .  Free  Volume  Additivi ty  Equations 

1 .  Introduction 

Of  the  simple  viscosity  theories  the  free  volume  treatment  has 
been  most  successful  in  accounting  for  the  temperature  dependence  of 
organic  liquid  viscosities  in  the  10”2  to  104  P  range.  According  to 
this  theory,  the  shear  viscosity  is  given  by  Eq.  (13): 

In  n  =  ln[(0.2/Bao2 )(MkT)1/2]  +  YV*/Vf  (13) 

The  free  volume  Vf  of  the  liquid  is  given  by 

Vf  =  VA«(T  -  To)  (39) 

where  V  is  the  specific  volume.  If  we  combine  Eqs.  (13)  and  (39), 
and  ignore  the  weak  temperature  dependence  due  to  the  T1/2  in  the  first 
term  on  the  right  hand  side  of  Eq.  (13),  we  obtain  an  expression  of  the 
form  of  the  VTF  Eq.,  (Eq.(35)). 

Granted  that  molecular  diameters,  ao ,  of  low  molecular  weight 
organic  molecules  are  all  of  comparable  magnitude,  Eq.  (13)  shows  that 
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Hence  we  will  introduce  the  further  simplification  that  all  of  the  B-j 
terms  in  Eqs.  (41)  and  (42)  are  the  same,  so  that 

In  n  =  EX-jAi  +  B/(T  -  ZXiTo-j)  (43) 

and 

1/n  =  LlXiexp(-Ai)]exp[-B/(T  -  EX^oi)] 

=  (zX^ai )exp[-B/(T  -  nXiTo-j )]  (44) 

where  ai  =  exp(-Ai)  in  Eq.  (44).  We  will  call  these,  the  master  viscosity 
Eqs.  They  can  be  justified  roughly  as  follows. 

We  have  mentioned  that  the  B  parameters  of  the  best  fits  to  the 
VTF  Eq.  in  Table  III  for  the  three  binaries  and  for  the  mixtures  of  all 
C1Q  hydrocarbons  are  nearly  the  same.  The  mean  values  of  B  are: 

581.26  for  ND,  NBC,  CD,  NBB,  DM0  and  all  their  mixtures 

536.26  for  NBC-DMO  mixtures 

637.48  for  CD-DMO  mixtures 

832.45  for  HNN-XTHDCPD  mixtures 

Using  these  B  values  fits  to  the  viscosities  of  the  respective  liquids 
could  be  obtained  which  were  only  slightly  worse  than  the  best  fits  of 
Table  III.  Parameter  for  VTF  Eq.  fits  using  the  mean  values  of  B  are 
given  in  Table  IX. 

Inspection  of  these  parameters  shows  that  forcing  a  constant 
value  of  B  orders  the  A  and  To  parameters,  such  that  the  A  and  To 
values  for  the  mixtures  are  intermediate  between  the  respective  A  and 
To  parameters  for  the  pure  components.  This  means  that  the  A  and  To 


parameters  for  any  of  the  liquids  are  given  approximately  by 

A  =  zXiA-j 


(45) 
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To  -  £XjTo-j  (45) 

where  subscript  i  refers  as  before  to  component  i. 

A,  and  To^  values  were  obtained  by  linear  least  squares  fits 
to  Eqs.  (45)  and  (46)  of  the  A  and  To  values  in  Table  IX.  For  the  five 
Cj0  hydrocarbons  these  fits  gave: 

AND  3  “7.701  To  nD  =  96.6 

aNBC  =  “7.505  To  NBC  =  107.8 

aCD  -  -7.003  To  qo  =  124.5  (47a) 

ADMO  =  -7.782  To  DM0  =97.5 

AnbB  3  -7.618  To  NBB  =  97.6 

STD.  DEV.  =  0.030  STD.  DEV.  =  1.9 

For  the  HNN-XTHDCPD  binary  the  fits  to  Eqs.  (45)  and  (46)  gave: 

AHNN  3  -7.077  Tohnn  =  152.3 

AXTHDCPD  =  -7-599  ToXTHDCPD  3  92.1  (47b) 

STD.  DEV.  =  0.015  STD.  DEV.  =  0.8 

Substituting  Eqs.  (45)  and  (46)  into  Eq.  (35)  justifies  the  serni- 
empi rical  Eq.  (43)  which  accounts  for  the  combined  composition/ 
temperature  dependence  of  the  viscosity. 

2 •  Fits  to  Master  Vlscosi ty  Equation  (4 3} 

We  have  carried  out  a  fit  to  Eq.  (43)  of  the  experimental 
viscosity  data  for  all  temperatures  and  compositions  of  each  of  the 
three  binary  systems  (NBC-DMO,  CD-DMO  and  HNN-XTHDCPD)  and  of  the 
combined  data  (1,  2,  3.  4,  and  5  -  component  liquids)  for  all  five  C 


hydrocarbons.  Eq.  (43)  is  linear  in  the  Ai  and  B  parameters,  so  that 
these  were  evaluated  by  a  linear  least  squares  computer  fit.  The  Toi 
parameters  were  set  equal  to  those  obtained  in  the  previous  section 
(Eq.  (47)).  The  relative  precision  or  %  error  for  all  viscosity 
determinations  are  about  the  same  and  Independent  of  the  magnitude  of 
n,  that  Is,  An/n  -  Ain  n  *  constant.  Hence  in  fits  to  Eq.  (43)  all 
values  of  In  n  were  weighted  equally.  The  parameters  for  the  fits  to 
Eq.  (43)  are  given  In  Table  X. 

In  Table  XI  a  comparison  Is  made  of  experimental  n  values  with 
those  calculated  via  Eq.  (43)  from  the  parameters  of  Table  X.  The 
quantity  tabulated  Is  100(ln  ncaic  -  In  nexp)  *h1ch  Is  approximately 
the  %  deviation  of  the  experimental  value  from  the  best  fit  value. 

Clearly  Eq.  (43)  does  not  fit  the  data  to  within  experimental  error. 

Pure  component  viscosities  are  underestimated  for  the  five  C1() 
hydrocarbons  and  overestimated  for  HNN  and  XTHDCPD.  Mixture  viscosities 
are  generally  overestimated  for  five  C^'s  and  underestimated  for  HNN- 
XTHDCPD.  The  disagreement  between  experimental  and  calculated  viscosities 
Is  largest  at  low  temperatures  and  where  the  differences  in  pure  component 
viscosities  are  largest.  That  Is,  the  fit  Is  better  for  the  NBC-DMO 
binary,  where  the  pure  component  viscosities  are  not  very  different, 
than  for  the  CD-DMO  and  HNN-XTHDCPD  binaries,  where  the  opposite  Is 
true. 

On  the  other  hand  Eq.  (43)  gives  a  remarkably  good  account  of 
the  data.  For  systems  of  N  components  it  describes  the  viscosities  of 
both  pure  liquids  and  mixtures  as  a  function  of  temperature  using  (N  -  1) 
fewer  adjustable  parameters  than  are  required  by  the  VTF  Eq.  to  account 
for  n  as  a  function  of  T  of  the  pure  liquids  alone.  For  the  HNN-XTHDCPD 
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binary  system,  for  example,  the  experimental  data  cover  a  range  in  In  n 
of  8.61.  The  standard  deviation,  in  In  n  from  the  fit  to  Eq.  (43)  is 
0.056,  i.e.,  only  0.7%  of  the  range  in  In  n.  The  error  in  In  n  for  the 
worst  data  point  is  0.15,  only  1.7%  of  the  range  in  In  n. 

3.  fits  to  Master  Viscosity  Equation  (44) 

Fits  of  the  experimental  viscosity  data,  similar  to  those 
described  in  the  previous  section,  were  also  carried  out  for  Eq.  (44). 
The  Toi  values  were  again  set  equal  to  those  of  Eq.  (47).  Eq.  (44)  is 
not  linear  in  B,  so  the  best  fit  B  parameter  was  obtained  by  trial  and 
error.  For  each  trial  value  of  B  the  corresponding  a^  values  were 
obtained  from  a  linear  least  squares  fit.  As  mentioned  in  the  previous 
section,  the  experimental  error  in  each  value  of  In  n  is  approximately 
the  same.  Consequently,  the  error  in  the  fluidity,  &(l/n),  is  not 
constant  but  is  proportional  to  1/n.  Hence  in  carrying  out  the  fits  to 
Eq.  (44)  each  fluidity  point  was  weighted  by  the  fluidity  Itself.  That 
is,  the  ai  and  B  parameters  were  chosen  to  minimize  the  following 
quantity: 

x2*  ‘  ^j.calc)2^1^2 

where  n  is  the  number  of  data  points  Included  in  the  fit. 

The  parameters  for  the  best  fits  to  Eq.  (44)  are  also  given  in 
Table  X.  In  Table  XII  are  tabulated  values  of  100(ln  neXp  -  In  ncaic) 
for  Eq.  (44). 

For  the  CD-DMO  binary  system  Eq.  (44)  gives  a  markedly  better 
fit  to  the  data  than  does  Eq.  (43).  For  the  other  two  binaries  and 
for  the  combined  data  for  all  five  C10  hydrocarbons  there  is  not  a 
great  deal  of  difference  between  the  goodness  of  the  fits  obtained  from 
the  two  master  viscosity  equations. 
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4.  Free  Volume  Interpretation  of  Viscosity 
The  free  volume  theory  assumes  that  a  molecule  can  undergo 
diffusive  transport  only  when  the  free  volume  in  its  vicinity  is  equal 
to  or  greater  than  a  certain  critical  value  v*.  The  shear  viscosity  is 
given  by  the  Eq.  (14): 

In  n  -  ln[(l/eao2 Mkm/Sir)1/2]  ♦  (l/2)ln  T  +  {yv*/[va«(T  -  To)]} 

=  A'  +  (l/2)ln  T  +  B'/(T  -  T  ) 

If  we  equate  this  to  the  VTF  Eq.,  In  n  =  A  +  B/(T  -  To)  and  also  equate 
the  temperature  derivations  of  the  two  expressions 

din  n/dT  ■  -B/(T  -  T2)2  =  1/2T  -  B'/(T  -  T2)2 

we  get: 

A  -  ln[(1/8ao2  )(km/8ir)1/2]  +  (1/2)[1  +  In  T  -  (T2/T)]  (48) 

B  -  yvWAo  -  (T  -  T2)2/2T  (49) 

Over  the  temperature  range  covered  by  this  study  the  second  term  on  the 
right  hand  side  of  Eq.  (48)  Is  a  very  weak  function  of  temperature, 
while  the  second  term  on  the  right  hand  side  of  Eq.  (49)  Is  quite  small 
compared  to  our  experimental  B  values.  For  example  1/2(1  +  In  T  -  T2/T) 
for  T2  ~w100K  changes  from  2.98  to  3.23  between  -50  and  50°C,  compared 
to  a  typical  value  of  A  of  arou,..'  -7  (Table  X).  (T  -  T2)2/2T  for  T2  -• 
100K  Is  about  77  for  T  ■  50°C  and  about  34  for  T  ■  -50°C,  compared  to 
B  -  600. 

The  molecular  diameters  aoi  calculated  at  0°C  for  five  CJo 
hydrocarbons,  HNN  and  XTHDCPD  from  Eq.  (48),  and  the  A}  parameters  of 
Eq.  (47)  using  0.2  for  the  value  of  8  are  shown  in  Table  XIII. 
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These  are  compared  with  and  seen  to  be  of  the  same  magnitude  as  the  mean 
molecular  diameters  calculated  from  the  densities  at  0°C,  (M/p0Na)5/3, 
where  M  is  molecular  weight  and  Na  Avogadro's  number. 

These  molecules  are  not  spherical.  To  obtain  alternative 
estimates  of  their  dimensions  we  constructed  scale  frame  work  molecular 
models  of  them.  The  molecular  model  set  is  manufactured  by  Prentice  - 
Hall,  Inc.,  Englewood  Cliffs,  N.J.  By  examining  the  projections  of 
each  model  on  a  background  area  when  viewed  from  different  aspects,  we 
determined  the  dimensions  (L,  W,  H)  of  a  right  parallopiped  which  would 
just  enclose  the  entire  molecule;  these  are  also  listed  in  Table  XII. 

The  CD,  HNN  and  XTHDCPO  molecules  are  rigid,  but  ND,  NBC,  DM0  and  NBB 
are  flexible  and  passes  a  number  of  configurations.  For  ND  the  molecular 
dimensions  were  averaged  over  four  different  configurations  expressed 
in  terms  of  the  staggered  (s)  or  eclipsed  (e)  arrangement  of  groups 
bonded  to  adjacent  carbons  along  the  backbone  of  the  molecule.  These 
four  ND  configurations  are  (1)  all  (s),  (2)  (e)  at  carbons  5-6, 

(3)  (e)  at  carbons  3  -  4  and  7  -  8  and  (4)  (e)  at  carbons  2  -  3,  5  -  6, 
and  8-9.  For  NBC  the  configuration  Is  such  that  the  n-butyl  group 
is  occupying  an  equatorial  position  on  the  cyclohexane  ring  and  the 
terminal -CH2CH3  group  and  cyclohexane  ring  are  In  the  trans  positions 
relative  to  each  other.  For  this  general  NBC  configuration  molecular 
dimensions  were  averaged  for  the  (s)  and  (e)  configurations  of 
groups  bonded  to  carbons  2  -  3  of  the  n-butyl  side  group.  For  DM0 
three  different  configurations  have  been  considered,  which  are  (1)  all 
staggered,  (2)  eclipsed  at  cartons  2  -  3,  6  -  7,  and  (3)  eclipsed  at 
carbons  1  -  2,  7  -  8.  Finally  for  NBB,  two  configurations  were  considered: 
(1)  all  staggered  and  (2)  eclipsed  at  carbons  2  -  3  of  n-butyl  group. 
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The  molecular  diameters  from  viscosity  a  j  are  remarkably  close  to 
realistic  molecular  diameters  from  the  models. 

They  correlate  best  with  the  maximum  molecular  dimensions  L 
obtained  from  the  models  as  shown  by  the  near  constancy  of  the  ao-j/L 
ratio. 

As  noted  in  Ref.  7,  the  constancy  of  the  aoi/L  values 
implies,  via  the  free  volume  theory,  that  the  critical  sized  void  into 
which  a  molecule  makes  a  diffusive  step  must  be  of  sufficient  size  to 
accomodate  the  longest  dimension  of  the  molecule.  This  suggests  in 
turn  that  the  molecules  are  freely  rotating  in  the  viscosity  range 
studied  here. 

The  differences  in  viscosities  of  the  C,„  molecules  studied  in 

1  o 

this  work  and  that  of  the  HNN  and  XTHDCPD,  are  mostly  due  to  differences 
both  in  the  Ai  and  To^  parameters.  From  Eqs.  (43)  and  (44),  what  is 
not  accounted  for  by  Ai  or  ai,  is  accounted  for  by  To*.  Of  the  five 
C1Q  liquids,  CD  and  DM0  differ  most  In  viscosity,  by  a  factor  of  3.28 
at  the  upper  temperature  limit  of  50°C  and  by  a  factor  of  10.79  at  the 
lower  limit  of  about  -55°C.  The  differences  In  the  Ai  of  CD  and  DM0 
(and  hence  the  differences  In  L )  account  for  a  factor  of  about  2.0 
difference  In  viscosity,  l.e.,  for  a  substantial  fraction  of  the  viscosity 
difference  in  the  experimental  range.  The  difference  between  the  Ai 
of  HNN  and  of  XTHDCPD  accounts  for  a  factor  of  about  1.6  difference  In 

i 

viscosities  of  these  components. 

In  view  of  the  simplicity  of  the  free  volume  model  for  viscosity 
of  liquids,  we  feel  that  this  Is  a  very  good  correlation  between 
molecular  structure  and  viscosity,  and  It  appears  then  that  the  free 
volume  model  gives  a  surprisingly  good  account  of  the  shear  viscosity  of 
these  liquids  and  their  mixtures  In  the  low  viscosity  range. 


.  Experimental  viscosities  and  densities  of  n-Decane 


Run  1 


T°(C) 

n  (P) 

p  (g/Cm3) 

49.90 

0.00611 

0.7075 

39.79 

0.00694 

0.7151 

30.06 

0.00790 

0.7223 

14.86 

0.00994 

0.7336 

0.00 

0.01291 

0.7447 

-  9.92 

0.01577 

0.7522 

-21.02 

0.02038 

0.7605 

-30.53 

0.02610 

0.7672 

-32.00 

(Frozen)-- 

-- 

Run  2 

T°(C) 

n  (P) 

P  (g/Cm3) 

49.88 

0.00613 

0.7075 

39.80 

0.00703 

0.7149 

30.50 

0.00788 

0.7218 

14.88 

0.00998 

0.7335 

-  0.08 

0.01294 

0.7450 

-10.24 

0.01585 

0.7518 

-20.33 

0.01996 

0.7594 

-31.31 

0.02670 

0.7676 
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Table  I.  Experimental  viscosities  and  densities  of  n-Butyl cyclohexane 


Run  1 


1  T°(C) 

n  (P) 

P  (g/Cm3 

|  50.24 

0.00829 

0.7767 

30.32 

0.01094 

0.7310 

1  14.74 

0.01424 

0.8028 

-  -  0.33 

0.01910 

0.8139 

*  -15.15 

0.02707 

0.8246 

|  -30.22 

0.04185 

0.8359 

-40.13 

0.05923 

0.8430 

{  -50.16 

0.08982 

0.8503 

-56.61 

0.12374 

0.8553 

i 

Run  2 

1 

|  T°(C) 

n  (P) 

p  (g/Cm3 

50.04 

0.00832 

0.7772 

|  30.25 

0.01102 

0.7912 

15.22 

t 

0.01414 

0.8026 

1  -  0.34 

0.01910 

0.8144 

•  -15.19 

0.02702 

0.8251 

'  -30.10 

0.04161 

0.8361 

[  -40.22 

1 

0.05936 

0.8436 

-50.63 

0.09213 

0.8512 

1  -58.51 

0.13670 

0.8570 
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Table  I .  Experimental  viscosities  and  densities  of  Ci s-Decahydro-naphthalene . 


Run  1 


1 

Jl(C) 

n_(P) 

p(g/Cm3) 

49.96 

0.01846 

0.8734 

1 

30.05 

0.02673 

0.8885 

1 

12.98 

0.03926 

0  .9017 

1 

-  0.09 

0.05596 

0.9118 

1 

-15.47 

0.08969 

0.9229 

f 

1 

-30.31 

0.15828 

0.9343 

-40.18 

0.24844 

0  .9418 

-50.25 

0.42821 

0.9494 

1 

1 

-54.77 

0.56464 

0.9533 

1 

s 

l 

Run  2 

! 

I 

T°(C) 

n  (P) 

p  (g/Cm3) 

49.53 

0.01852 

0.8726 

| 

30.28 

0.02660 

0.8865 

i 

-  0.77 

0.05687 

0.9100 

( 

-15.13 

0.08849 

0.9209 

i 

i 

-32.68 

0.17438 

0.9339 

-40.43 

0.25020 

0.9401 

< 

\ 

-50.57 

0.43170 

0.9482 

! 

-55.78 

0.59974 

0.9516 

i 
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Table  I 


.  Experimental  viscosities  and  densities  of  Dimethyloctane 


T°(C) 

50.03 

» 

39.97 

30.14 

14.98 

-  0.29 

-14.76 

-30.09 

-40.00 

-50.13 


nJPl 

0.00559 

0.00632 

0.00720 

0.00901 

0.01172 

0.01567 

0.02278 

0.03051 

0.04367 


p  (q/Cm3) 
0.7012 

0.7085 

0.7160 

0.7273 

0.7386 

0.7495 

0.7607 

0.7677 

0.7752 


-55.34 


0.05396 


0.7790 
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Tabl^e  I .  Experimental 

1 

viscosities  and  densities 

of  n-Butyl benzene 

\ 

i 

Run  1 

•  T°(C) 

JL.CP1 

p  (g/Cm3) 

49.92 

0.00694 

0.8360 

|  28.94 

0.00919 

0.8525 

-  0.14 

0.01470 

0.8753 

1  -15.65 

0.02044 

0.8876 

-32.38 

0.03238 

0.9010 

-40.32 

0.04202 

0.9070 

;  -50.31 

0.06192 

0.9148 

-53.61 

0.07173 

0.9173 

i 

1 

Run  2 

j  T°(C) 

n  <P) 

P  (g/cm3) 

49.90 

0.00696 

0.8366 

|  40.14 

0.0C786 

0.8441 

30.12 

0.00897 

0.8518 

'  15.11 

0.01126 

0.8635 

|  -  0.08 

0.01474 

0.8751 

-15.11 

0.02031 

0.8870 

{  -29.95 

0.03131 

0.8984 

-40.17 

1 

0.04332 

0.9070 

1  -50.18 

0.06205 

0.9144 

i  -54.84 

0.07652 

0.9179 

Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Decane  +  50%  n-Butylcyclohexane 


Run  1 


T°  (C ) 

n  (P' 

p  (g/Cm3 ) 

49.92 

0.00702 

0.7414 

39.84 

0.00799 

0.7485 

30.26 

0.00914 

0.7557 

15.39 

0.01152 

0.7666 

-  0.59 

0.01542 

0.7782 

-10.75 

0.01913 

0.7860 

-20.81 

0.02438 

0.7937 

-30.51 

0.03190 

0.8007 

-40.20 

0.04352 

0.8072 

-42.00 

Frozen  - 

-- 

Run  2 

T°(C) 

n  (P) 

P  (g/Cm3) 

49.91 

0.00705 

0.7407 

39.81 

0.00799 

0.7477 

29.65 

0.00920 

0.7553 

14.77 

0.01163 

0.7662 

-  0.03 

0.01524 

0.7771 

-10.52 

0.01897 

0.7849 

-20.69 

0.02453 

0.7925 

-30.41 

0.03267 

0.7999 

-40.35 

0.04474 

0.8067 

Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Decane  +  50%  cis-Decahydronaphthalene 


T°(C) 

nJPl 

iOaZCm3! 

49.86 

0.00927 

0.7829 

39.77 

0.01066 

0,7901 

29.96 

0.01242 

0.7974 

14.40 

0.01631 

0.8090 

-  0.27 

0.02188 

0.8198 

-10.73 

0.02788 

0.8277 

-20.32 

0.03584 

0.8346 

-30.33 

0.04852 

0.8420 

-40.19 

0.06835 

0.8488 

I 

t 


t 

.  Table  I 

l 

I 

I 
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Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Decane  +  50%  Oimethyloctane 


T°(C) 

n  (P) 

P  (g/Cm 

50.10 

0.00584 

0.7048 

39.97 

0.00662 

0.7121 

30.03 

0.00755 

0.7197 

15.06 

0.00944 

0.7307 

-  0,26 

0.01233 

0.7424 

-15.00 

0.01666 

0.7529 

-30.11 

0.02426 

0.7638 

-39.90 

0.03263 

0.7711 

-50.04 

0.04686 

0.7783 

40 


Table 


I. 


SKisss  :'5or:^yc„su,es  of  the  "•«-  °f 


T°(C) 

n  (P) 

P  (g/Cm3) 

50.03 

0.00622 

0.7642 

40.60 

0.00695 

0.7709 

30.67 

0.00792 

0.7790 

14.77 

0.00998 

0.7910 

-  0.56 

0.01293 

0.8026 

-15.07 

0.01730 

0.8135 

-29.96 

0.02481 

0.8248 

-40.10 

0.03345 

0.8322 

-44.66 

0.03899 

0.8359 
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Table  I .  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Butylcyclohexane  +  50%  cis-Decahydronaphthalene. 


T°(C) 

P  (9/Cm3) 

49.86 

0.01162 

0.8227 

30.00 

0.01 593 

0  .8379 

14.70 

0.02129 

0.8489 

-  0.37 

0.02974 

0.8599 

-14.73 

0.04338 

0.8708 

-30.36 

0.07177 

0.8824 

-41.77 

0.10861 

0.8902 

-51.37 

0.17644 

0.8978 

-55.18 


0.21745 


0.9007 


Table  I.  Experimental  viscosities  and  densities  of  the  mixture 
75%  n-Butylcyclohexane  +  25%  Dimethyl  octane 
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T°  (C) 

50.09 

40.21 

29.93 

15.04 

0.00 

-15.10 

-30.03 


-39.97 


-50.00 


-53.63 


n  (P) 

0.00739 

0.00848 

0.00982 

0.01248 

0.01652 

0.02317 

0.03501 

0.04907 

0.07336 

0.08641 


p  (9/Cm3) 

0.7506 

0.7637 

0.7706 

0.7824 

0.7929 

0.8042 

0.8149 

0.8223 

0.8290 

0.8317 
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Table_I -  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Butyl cyclohexane  +  50%  Dimethyl  octane 


T°_(C) 

n  (P) 

p  (g/Cm3) 

50.21 

0.00673 

0.7379 

40.16 

0.00766 

0.7452 

30.16 

0.00880 

0.7523 

15.26 

0.01111 

0.7641 

-  0.09 

0.01466 

0.7748 

-14.84 

0.02014 

0.7853 

-29.97 

0.03008 

0.7963 

-39.89 

0.04132 

0.8030 

-50.00 

0.06073 

0.8102 

-55.96 

0.07950 

0.8146 
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Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
25%  ri-Buty  lcyc  lohexane  +  75%  Dimethyloctane 


T°(C) 

50.07 

40.19 

30.07 

14.89 

0.04 

-15.22 

-30.10 

-40.01 

-50.56 


n  (P) 

0.00613 

0.00694 

0.00794 

0.00994 

0.01295 

0.01782 

0.02603 

0.03545 

0.05232 


P  (g/Cm3 ) 

0.7194 

0.7266 

0.7339 

0.7455 

0.7564 

0.7675 

0.7778 

0.7857 

0.7934 


-54.65 


0.06238 


0.7960 
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Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  n-Butyl benzene  +  50%  n-Butylcyclohexane 


jajn 

JL  (flACm? ) 

49.98 

0.00726 

0.8040 

30.06 

0.00944 

0.8196 

14.84 

0.01192 

0.8307 

-  0.41 

0.01575 

0.8430 

-15.61 

0.02175 

0.8544 

-30.69 

0.03239 

0.8658 

-40.68 

0.04443 

0.8728 

-50.55 

0.06457 

0.8806 

-52.45 


0.06999 


0.8822 


Table  1.  Experimental  viscosities  and  densities  of  the  mixture  of 
25',  i  Dimethyl  octane  +  75vi  cis-Decahydronaphthalene 


^(0 

50.03 

40.01 

30.01 

14.75 

-  0.22 

-15.13 

-29.96 

-40.21 

-50.10 


n  (P) 

0.01239 

0.01450 

0.01722 

0.02311 

0.03254 

0.04877 

0.08003 

0.12045 

0.19201 


o  (g/Cm3 ) 

0.8244 

0.8314 

0.8390 

0.8501 

0.8616 

0.8725 

0.8836 

0.8912 

0.8984 


-54.38 


0.24133 


0.9013 
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Table  1-  Experimental  viscosities  and  densities  of  the  mixture  of 
50%  ci s-Decahydronaphthal ene  +  50t  Dimethyl  octane 


T°(C) 

n(P) 

p  (g/Cm3) 

50.04 

0.00906 

0.7796 

40.25 

0.01044 

0.7867 

30.10 

0.01211 

0.7939 

14.82 

0.01 585 

0.8050 

-  0.10 

0.02147 

0.8161 

-15.20 

0.03097 

0.8271 

-30.15 

0.04811 

0.8380 

-40.00 

0.06854 

0.8447 

-50.39 

0.10550 

0.8523 

-55.26 


0.13530 


0.8557 


Table 


1 


Experimental  viscosities  and  densities  of  mixtures  of 
75%  Dimethyloctane  +  25%  cis-Decahydronaphthalene 


I 


l 

I 

I 

I 

I 

I 

I 

I 

I 

! 

I 

I 

( 

! 


T°(C) 

50.11 

40.12 

29.87 

15.28 

-  0.38 

-15.02 


-30.16 


-39.91 


-50.18 


-53.33 


n  (P)  n  (9/Cni3) 


0.00696 

0.00855 

0.00918 

0.01159 

0.01545 

0.02133 

0.03206 

0.04405 

0.06528 

0.07481 


0.7386 

0.7460 

0.7536 

0.7644 

0.7759 

0.7865 

0.7974 

0.8052 

0.8122 

0.8144 
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Table  I.  Experimental  viscosities  arid  densities  of  the  mixture  of 
50%  ci s-Oecahydronaphthalene  +  50%  n-Butylbenzene 


t!Ic) 

50.14 

29.98 

14.51 

-  0.12 

-14.87 

-30.79 

-40.25 

-50.25 


0.01000 

0.01309 

0.01699 

0.02268 

0.03229 

0.05118 

0.07159 

0.10889 


p-la/C'i 

0.8534 

0.8692 

0.8810 

0.8920 

0.9037 

0.9157 

0.9226 

0.9305 


-53.72 


0.12891 


0.9334 


BO 


Table  1.  Experimental  viscosities  and  densities  of  the  misture  of 
501  i  n-Butyl benzene  +  50%  Dimethyl  octane 


tec) 

50.24 

40.11 

30.06 

15.11 

0.05 

-14.94 

-30.00 

-40.04 

-50.52 


-54.68 


.nJ.P) 

0.00597 

0.00673 

0.00766 

0.00949 

0.01222 

0.01640 

0.02360 

0.03175 

0.04572 

0.05387 


oJq/Cm3J 

0.7608 

0.7688 

0.7762 

0.7875 

0.7990 

0.8106 

0.8218 

0.8296 

0.8375 

0.8409 


Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
33.3%  n-Decane  +  33.3%  n-Butycyclohexane  +  33.3%  cis- 
Decahydronaphthalene 

T°(C)  ilPj  p  (g/Cm3) 

49.91  0.00892  0.7811 

30.26  0.01187  0.7951 

15.24  0.01531  0.8063 

-  0.23  0.02085  0.8177 

-15.11  0.02963  0.8288 

-30.65  0.04642  0.8402 

-40.29  0.06503  0.8473 


-50.41 


0.09862 


0.8545 
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Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 
33.31'  n-Buty  1  benzene  +  33. 3%  n-Butylcyclohexane  +  33.3'*. 
n-Decane 


T'^(C) 

49.72 

39.75 

30.11 

14.79 

-  0.17 

-14.85 

-29.93 

-40.01 

-50.52 


i)  (P) 

0.00679 

0.00768 

0.00873 

0.01103 

0.01433 

0.01941 

0.02851 

0.03885 

0.05619 


P_(g/Cm3) 

0.7686 

0.7763 

0.7831 

0.7951 

0.8064 

0.8169 

0.8285 

0.8359 

0.8438 


-54.56 


0.06653 


0.8465 


5.3 


Tablet.  Experimental  viscosities  and  densities  of  the  mexture  of 

33. 3^  ci s-Decahydronaphthalene  +  33.3%  n-Butyl benzene  *  33.3% 
n-Decane 


Illc) 

_nJP_) 

P  (g/Cm3) 

49.73 

0.00806 

0.7979 

39.87 

0.00918 

0.8056 

30.10 

0.01054 

0.8126 

14.88 

0.01339 

0.8240 

-  0.09 

0.01771 

0.8354 

-15.35 

0.02475 

0.8466 

-32.88 

0.04021 

0.8601 

-40.31 

0.05119 

0.8651 

-50.54 


0.07564 


0.8729 
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Table  I.  Experimental  viscosities  and  densities  of  the  mixture  of 

33.3"  n-Butylbenzene  +  33.3%  n-Butylcyclohexane  +  33.3%  cis- 
Decahydronaphthal ene 


fjc) 

n  _  ( P ) 

p  (g/Cnr 

50.04 

0.00917 

0.8265 

30.19 

0.01218 

0.8414 

14.87 

0.01573 

0.8526 

0.04 

0.02106 

0.8644 

-15.14 

0.03001 

0.8753 

-30.13 

0.04594 

0.8870 

-40.30 

0.06538 

0.8945 

-49.97 

0.09736 

0.9018 

-51.07 

0.10242 

0.9026 

-54.94 


0.12284 


0.9051 
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Table  l-  Experimental  viscosities  and  densities  of  the  mixture  of 
25%  n-Decane  +  25%  n-Butyl cyclohexane  +  25%  ci s-Decahydro- 
naphthalene  +  25%  n-Butyl benzene 


T°(C) 

n  (P) 

p  (g/Cm3) 

49.93 

0.00809 

0.7925 

39.86 

0.00926 

0.8001 

30.13 

0.01062 

0.8071 

15.41 

0.01345 

0.8181 

-  0.07 

0.01791 

0.8297 

-14.92 

0.02492 

0.8406 

-30.15 

0.03765 

0.8521 

-40.03 

0.05200 

0.8594 

-50.03 

0.07645 

0.8671 

-54.87 

0.09449 

0.8710 
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Tab  I  e 


Experimental  viscosities  and 
20«  of  5  C  components 


Run  1 

T°(c) 

n  (P) 

50.05 

0.00747 

40.24 

0.00848 

30.17 

0.00975 

14.88 

0.01242 

-  0.10 

0.01635 

-14.91 

0.02257 

-30.00 

0,03373 

-39.89 

0.04639 

-49.88 

0.06786 

-55.66 

0.08762 

Run  2 

T°(C) 

n  (P) 

51.83 

0.00731 

41.56 

0.00833 

30.91 

0.00965 

9.34 

0.01369 

-  1.68 

0.01686 

-12.68 

0.02141 

-24.68 

0.02897 

-37.12 

0.04227 

-49.20 

0.06588 

-54.44 

0.08328 

densities  of  the  mixture  of 


p  (g/Cm3) 
0.7725 
0.7794 
0.7868 
0.7985 
0.8096 
0.8206 
0.8319 
0.8390 
0.8463 
0.8504 

P  (g/Cm3 ) 
0.7714 
0.7790 
0.7868 
0.8029 
0.8113 
0.8197 
0.8284 
0.8372 
0.8463 
0.8503 


! 


Table  I.  Experimental  viscosities  and  densities  of  HNN 


I 

I 

I 

I 

I 

I 

I 

! 

I 

I 

I 

( 

( 

I 


T°(C) 

49.55 

39.74 

29.85 

14.94 

0.04 

-10.26 

-15.12 

-21.69 

-29.79 

-39.93 

-50.50 


n  (P) 

0.10568 

0.14357 

0.20099 

0.36634 

0.76796 

1.44592 

2.02257 

3.38706 

7.02662 

21.73532 

96.66266 


p  (g/cm3 ) 

1.0690 

1.0760 

1.0831 

1.0941 

1.1048 

1.1128 

1.1169 

1.1215 

1.1274 

1.1358 

1.1443 


-53.61 


Frozen 


Table  I.  Experimental  viscosities  and  densities  of  mixture  of 
80%  HNN  +  20%  XTHDCPD 


T°(C) 

n  (P) 

p  (g/Cm3 ) 

50.00 

0.07283 

1.0415 

29.90 

0.12759 

1.0551 

15.28 

0.21149 

1.0664 

-  0.04 

0.40261 

1.0791 

-14.94 

0.88599 

1.0903 

-30.07 

2.53050 

1.1018 

-40.09 

5.77212 

1.1094 

-50.39 

18.32347 

1.1171 

-54.92 

33.46471 

1.1207 

Table 


Experimental  viscosities  and  densities  of 
60%  HNN  +  40 %  XTHDCPD 


T°(C) 

49.85 

30.09 

15.10 

-  0.04 

-14.74 

-30-08 

-39-99 

-50.52 


n  (P) 

0.05071 

0.08266 

0.12848 

0.22233 

0.42766 

0.98465 

2.01006 

5.00710 


ixture  of 

P  (g/Cma ) 

1.0132 

1.0281 

1.0396 

1.0509 

1.0618 

1.0730 

1.0810 

1.0890 


-54.00 


7.08425 


1.0926 
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Table  I.  Experimental  viscosities  and  densities  of  mixture  of 
40%  HNN  +  60%  XTHDCPD 


T°(C) 

n  (P) 

p  (g/Cm3) 

49.92 

0.03525 

0.9821 

30.03 

0.05451 

0.9965 

14.62 

0.08201 

1.0087 

-  0.02 

0.12821 

1.0196 

-15.00 

0.22525 

1.0313 

-30.06 

0.45584 

1.0433 

-41.40 

0.87379 

1.0522 

-50.52 

1.64076 

1.0584 

-53.73 

2.10577 

1.0607 
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1.  Experimental  viscosities  and  densities  of  mixture  of 
20%  HNN  +  80,:  XTHDCPD 


.H(C)_ 

50.00 

30.15 

14.88 

0.00 

-15.00 

-30.00 

-40.08 

-50.00 


n  (P) 

0.02490 

0.03662 

0.05203 

0.07778 

0.12451 

0.22295 

0.35547 

0.61178 


p  (9/Cnr) 

0.9492 

0.9640 

0.9759 

0.9870 

0.9983 

1.0099 

1.0175 

1.0251 


-51.51 


0.69513 


1.0262 
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Tab  I  e 


.!_>  Experimental  viscosities  and  densities  of 


T°jC) 

49.99 

30.05 

15.38 

0.00 

-14.88 

-30.07 


-39.60 


-50.20 


-54.19 


ri  (P) 

0.01765 

0.02488 

0.03369 

0.04813 

0.07401 

0.11712 

0.16979 

0.27131 

0.33199 


XTHDCPD 

p  (g/Cm3) 
0.9126 

0.9273 

0.9392 

0.9508 

0.9618 

0.9735 

0.9808 

0.9892 

0.9914 
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!  Table  I  la. 

i 

i 

|  Molar 

|  ND  (Run 

(Run 

I  NBC  (Run 

j  (Run 

CD  (Run 

|  ( Run 

DM0 

NBB  (Run 
NBB  (Run 
0.5ND  - 

|  0.5ND  - 

0.5ND  - 

I  0.5MD  - 

|  0.5NI)  - 

0.5HBC  - 

|  0.75NBC 

0.5NBC  - 

I  0.25NBC 

j  0.5NBC  - 

0.75CD  ■ 

I 


Parameters  tor  density  equations  for  ND,  NBC,  CD,  DM0, 
NBB,  and  their  mixtures 


p  (g/Cm1)  =  p(o)  -  bT(°C) 


Compos i ti on 

T  range(°C) 

p(o) 

1 0*  b 

Std. 

1) 

-31 

to 

50 

0.7447 

7.44 

0.0001 

2) 

-31 

to 

50 

0.7444 

7.41 

0.0002 

1) 

-57 

to 

50 

0.8135 

7.35 

0.0001 

2) 

-59 

to 

50 

0.8139 

7.37 

0.0002 

1) 

-55 

to  50 

0.9114 

7.59 

0.0002 

2) 

-56 

to 

50 

0.9096 

7.54 

0.0002 

-56 

to 

50 

0.7383 

7.39 

0.0002 

1) 

-55 

to 

50 

0.8752 

7.77 

0.0001 

2) 

-54 

to 

50 

0.8753 

7.87 

0.0001 

0.5NBC  (Run  1) 

-40 

to 

50 

0.7780 

7.36 

0.0002 

0.5NBC  (Run  2) 

-40 

to 

50 

0.7776 

7.38 

0.0001 

0.5CD 

-40 

to 

50 

0.8195 

7.36 

0.0002 

0.5DM0 

-50 

to 

50 

0.7417 

7.35 

0.0002 

0.5NBB 

-45 

to 

50 

0.8020 

7.58 

0.0002 

0.5CD 

-55 

to 

50 

0.8598 

7.39 

0.0002 

-  0.25DM0 

-54 

to 

50 

0.7922 

7.53 

0.0015 

0.5DM0 

-56 

to 

50 

0.7743 

7.23 

0.0003 

-  0.75DM0 

-54 

to 

50 

0.7562 

7.35 

0.0002 

•  0.5NBB 

-53 

to 

50 

0.8423 

7.60 

0.0003 

■  0.25DM0 

-54 

to 

50 

0.8613 

7.41 

0.0002 
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0.5CD  -  0.5DM0 

-55 

to 

50 

0.8158 

7.24 

0.0002 

0.25CD  -  O.75DM0 

-53 

to 

50 

0.7755 

7.33 

0.0002 

0.5CD  -  0.5NBB 

-54 

to 

50 

0.8920 

7.67 

0.0002 

0.5DM0  -  0.5NBB 

-55 

to 

50 

0.7991 

7.61 

0.0001 

1 / 3ND- 1 / 3NBC- 1 / 3C0 

-50 

to 

50 

0.8176 

7.35 

0.0001 

1/3ND-1/3NBC-1/3NBB 

1 /3ND-1 /3CD-1/3NBB 
1/3NBC-1/3CD-1/3NBB 

-55 

-51 

-55 

to 

to 

to 

50 

50 

50 

0.8059 

0.8352 

0.8641 

7.49 

7.47 

7.52 

0.0002 

0.0002 

0.0002 

0.25ND  -  0.25NBC  -  0.25CD 
-0.25NBB 

-55 

to 

50 

0.8297 

7.47 

0.0002 

0. 2ND  -  0.2NBC  -  0.2CD 
-0.2DM0  -  0.2NBB 

-56 

to 

50 

0 . 8(094 

7.41 

0.0002 

6 


Table  1 1 b .  Parameters  for  density  equations  for  HNN,  XTHDCPD  and 
thei r  mixtures 


p  (g/Cni3)  =  p(o)  -  bT(°C) 


Molar  Composition 

T  ranqe(°C) 

p(o) 

1 04  b 

Std. 
Dev,  p 

HNN 

-50  to  50 

1.1056 

7.50 

0.0005 

0.8HNN  -  0.2XTHDCPD 

-55  to  50 

1.0788 

7.62 

0.0005 

0.6HNN  -  0.4XTHDCPD 

-54  to  50 

1.05085 

7.56 

0.0004 

0.4HNN  -  0.6XTHDCPD 

-54  to  50 

1.0199 

7.65 

0.0004 

0.2HNN  -  0.8XTHDCPD 

-51  to  50 

0.9871 

7.60 

0.0001 

XTHDCPD 

-54  to  50 

0.9506 

7.61 

0.0003 

I 

j 

j 

0'.' 

iahle  Ilia.  Host,  fit  V T F  equation  parameters  lor  shear  viscosity  of 
NO,  NBC.  (.11.  DMO,  NUB,  anti  their  mixtures 


In  O’)  -  A  +  B/[T(K)  -  To] 

I 


|  Molar  Compos i  tiori 

T  range('C) 

A 

B 

T2_ 

Std. 
Dev.  In 

£  ill)  (Run  1) 

-31 

to 

50 

-7.7764 

614.11 

94 

0.001 

(Run  2) 

-31 

to 

50 

-7.7442 

605.03 

95 

0.002 

(  NBC  (Run  I) 

-57 

to 

50 

-7.3986 

549.44 

113 

0.004 

-  (Run  2) 

-59 

to 

50 

-7.3679 

541.14 

114 

0.U04 

"  CD  (Run  1) 

-55 

to 

50 

-7.6572 

796.29 

106 

0.004 

|  (Run  2) 

-56 

to 

50 

-7.6540 

795.18 

106 

0.005 

DMO 

-56 

to 

50 

-7.6158 

534.10 

104 

0.005 

|  MBB  (Run  1) 

-55 

to 

50 

-7.2973 

483.68 

116 

0.003 

(Run  2) 

-54 

to 

50 

-7.1709 

451.68 

120 

0.001 

■  0.5ND  -  0.5NBC  (Run  1) 

-40 

to 

50 

-7.6013 

584.93 

102 

0.002 

|  0.5ND  -  0.5NBC  (Run  2) 

-40 

to 

50 

-7.6843 

613.23 

98 

0.003 

0. 5N D  -  0.5CD 

-40 

to 

50 

-7.6322 

663.30 

98 

0.004 

|  0.5ND  -  0.5DM0 

-50 

to 

50 

-7.6243 

547.89 

103 

0.004 

0.5ND  -  0.5U8B 

-45 

to 

50 

-7.5610 

554.51 

100 

0.003 

1  0.5NBC  -  0.5CD 

-55 

to 

50 

-7.5108 

658.00 

108 

0.003 

|  0. 75NBC  -  0.2 5 DMO 

-54 

to 

50 

-7.4807 

550.85 

110 

0.005 

0.5NBC  -  0.5DM0 

-56 

to 

50 

-7.4730 

529.26 

no 

0.005 

|  0.25NBC  -  0.75DMO 

-55 

to 

50 

-7.5039 

517.64 

109 

0.004 

0.5HBC  -  0.5NBB 

-53 

to 

50 

-7.4491 

549.18 

106 

0.004 

1  0.7‘jCD  -  0.25DM0 

-54 

t.0 

50 

-7.5650 

686.71 

107 

0.002 

I 


$ 


0.5CD  -  0.5DM0 

-55 

to 

50 

-7.5067 

604.51 

108 

0.003 

0.25CD  -  0 . 7 SUMO 

-53 

to 

50 

-7.5645 

565.81 

106 

0.004 

0.5CD  -  0.5NBB 

-54 

to 

50 

-7.2042 

543.55 

114 

0.003 

0.5 DMO  -  0.5NBB 

-55 

to 

50 

-7.4655 

511.02 

106 

0.004 

1/3ND-1/3NBC-1/3CD 

-50 

to 

50 

-7.5518 

621.52 

104 

0.003 

1 / 3ND- 1 / 3NBC- 1  / 3NBB 

-55 

to 

50 

-7.5352 

562.32 

102 

0.003 

1/3ND-1/3CD-1/3NBB 

-51 

to 

50 

-7.4812 

586.03 

103 

0.002 

1/3NBC-1 /3CD-1/3MBB 

-55 

to 

50 

-7.4028 

584.63 

108 

0.004 

0.25ND  -  0.25NBC 
-  0.25CD  -  0.25NBB 

-55 

to 

50 

-7.4598 

577.65 

105 

0.003 

0. 2ND  -  0.2NBC  -  0.2CD 
-0.2DM0  -  0.2NBB 

-56 

to 

50 

-7.4569 

554.62 

107 

0.004 

I 


able  111b-  Bos t  fit  VII  equation  parameters  for  shear  viscosity  of 
MNN,  XrilijCPfj  ami  their  mixtures 


In  ,,((>)  -  A  t  B/LT(K)  -  T„] 


Mo  l  ar  Coiiij/os  i  t  i  on 

T  ramje('_'Cj  A 

B 

T, 

sr.d. 
Dev . 

HMN 

-50  to  50  -7.1247 

838.83 

151 

0.004 

0.8HNN  -  0.2XTIIDCPD 

-55  to  50  -7.1276 

821.79 

141 

0.005 

0.61  INN  -  0. 4XTHDCPD 

-54  to  50  -7.1661 

804.58 

131 

0.004 

0.4HNN  -  0.6XTHDCPD 

-54  to  50  -7.3421 

820.45 

118 

0.002 

0.2HNN  -  0.8XTHDCPD 

-50  to  50  -7.5266 

845.24 

103 

0.004 

XTHDCPD 

-54  to  50  -7.7096 

864.80 

88 

0.006 
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Table  IVa.  Comparison  of  densities  obtained  in  the  present  investigation 
with  those  reported  in  the  literature. 


p  (g/Cm3) 


Compound 

T(UC) 

Present 

Results 

Literature 

ND 

0 

0.744r, 

0.7444 

0. 74487a ,  0. 17448b 

20 

0.7298, 

0.7296 

0.729963,  0. 7298b  ,  0.73014C, 
0. 72994d ,  0.72987e 

30 

0.7224, 

0.7222 

0. 72238a 

NBC 

20 

0.7988, 

0.7992 

0.79918f 

30 

0.7915, 

0.7918 

0.791 76f 

CD 

0 

0.9114 

0.91209 

20 

0.8962 

0.89679 ,  0 . 8967 1 h 

30 

0.8886 

0.88929,  0.8891  lh 

dBruylants  (unp.).  Tinmermans. 
bChavanne  and  Took  (1932).  Timmermans. 
cShepard,  Henne  and  Midgley  (1931).  Timmermans. 
dBruun  and  Hicks  -  Brunn  (1932).  Timmermans. 
eMears,  Fooksonetal  (1950).  Timm.  Volume  II. 
^Forziati  and  Rossini  (1949).  Timm.  Volume  II. 
9Seyer  and  Davenport  (1941).  Tinmermans. 

hCamin  and  Rossini  (1955).  Timm.  Volume  II. 


7u 


MO 

0 

0.7383 

0.7 3785 1 

15 

0.7272 

0.726401 

30 

0.7161 

Q.714941 

NBB 

-45.2 

0.9103,  0.9103 

O.O^IO^ 

-22.9 

0.8930,  0.8933 

0.89440J 

0 

0.8752,  0.8753 

0.87608j,  0.87609"’ 

20 

0.8597,  0.8596 

0.86013f,  0.8601 3k ,  0.8603 

40 

0.8442,  0.8438 

0.84389k 

’Timmermans  and  Hennaut  -  Roland  (1929).  Timmermans. 

JMassart  (1936).  Timmermans 

^Donaldson  and  Quayle  (1950).  Timm.  Volume  II. 

^irch,  Deanetal  (1949). 

mB.  E.  (unp.).  Timmermans. 


Ttible  IVb.  Comparison  of  viscosities  obtained  in  present  investigation 
with  those  reported  in  the  literature. 


n  (CP) 

Cornpound 

IC'CJ 

Present _Resu 1 1  s 

Literature 

ND 

25 

0.850,  0.851 

0 . 852  7a 

NBC 

-55 

11.38,  11.39 

11. 7b 

-25 

3.57,  3.56 

3. 58b 

0 

1.891,  1.892 

1.90b 

25 

1.189,  1.189 

1 .205b 

50 

8.360,  8.356 

8 . 56 1 b 

CD 

0 

5.54,  5.52 

5.620c 

10 

4.23,  4.22 

4. 300C 

20 

3.33,  3.32 

3. 381 C 

30 

2.683,  2.677 

2.723C 

DM0 

15 

0.896 

0.889d 

30 

0.720 

0.703d 

Shepard,  Henne  and  Midgley  (1931).  Timmermans . 

^Carpenter  Davies  and  Matheson,  J.  Chem.  Phys.,  46,  2451  (1966). 

cSeyer  and  Leslie  (1942).  Timmermans. 

dTimmermans  and  Hennaut  -  Roland  (1929).  Timmermans. 
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MOB 


-50.7 

6.37,  6.32 

6 . 39e 

-30.5 

3.09,  3.06 

3. 12e 

-10.3 

1.825,  1.815 

1.82e 

9.8 

1.228,  1.229 

1 . 2  3e 

30.0 

0.898,  0.905 

0.897' 

Barlow,  Lamb  and  Matheson,  Proc.  Roy.  Soc.,  A292,  322  (1966) 
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Table  Va .  Test  of  isothermal  additivity  rules  for  shear  viscosities 
of  mixtures  of  Cl(J  hydrocarbons.  nex  is  the  experimental 
mixture  of  viscosity. 


Molar  Coiiifjosj  tiort 

(tXi/ni) 

50  'C 

1  /"ex 
-3jy;c 

exp(r.X-j  In 
50"  C 

ni )/"ex 
-_3_0l'C 

0.5ND  -  0.5NBC 

1.00 

1.01 

1.02 

1.04 

0.5ND  -  0.5CD 

1.00 

0.95 

1.15 

1.36 

0. 5ND  -  0.5DM0 

1.00 

0.99 

1.00 

1.00 

0.5ND  -  0.5NBB 

1.05 

1.11 

1.05 

1.12 

0.5NBC  -  0.5CD 

0.99 

0.93 

1.07 

1.14 

0.75MBC  -  0.2 5 DM0 

0.99 

0.98 

1.01 

1.02 

0.5NBC  -  0.5DM0 

0.99 

0.98 

1.01 

1.02 

0.25NBC  -  0.75DM0 

0.99 

0.99 

1.01 

1.02 

0.5NBC  -  0.5NBB 

1.04 

1.10 

1.05 

1.12 

0.75CD  -  0.25DM0 

0.95 

0.79 

1.11 

1.21 

0.5CD  -  0.5DM0 

0.95 

0.83 

1.12 

1.25 

0.25CD  -  0.75DM0 

0.97 

0.91 

1.08 

1.15 

0.5CD  -  0.5NBB 

1.02 

1 .02 

1.14 

1.38 

0.50M0  -0.5NBB 

1.04 

1.10 

1.04 

1.11 

1/3ND-1/3NBC- 1/3CD 

0.99 

0.95 

1.10 

1.21 

1/3ND-1/3NBC-1/3NBB 

1  .04 

1.09 

1.05 

1.11 

1 / 3ND- 1 / 3CD- 1 / 3NBB 

1 .04 

1.04 

1.15 

1.34 

1/3NBC-1/3CD-1/3HBB 

1.03 

1 .04 

1.12 

1.26 

0.25ND  -  0.25NBC 
-  0.25CD  -  0.25NBB 

1 .03 

1.04 

1.11 

1  .26 

0.2NU  -  0.2NBC  -  0.2C0 
-  0.2DM()  -  0.2NBB 

1.02 

1 .01 

1.10 

1.21 

Tab  Ip  _Vb.  Tests  of  isothermal  additivity  rules  for  shear  viscosities 
mixtures  of  HNN  and  XTHDCPD.  eex  is  the  experimental 
mixture  viscosity. 


( >;  X  ^  /  n  -j ! 

!  /''ex 

exp(>:Xi  In 

c i  )/r'ex 

Molar  Composition 

50  “  C 

-3Q°C 

50^C 

-30°  C 

0.8HNN  -  0.2XTHUCPI) 

0.72 

0.22 

1.01 

1.27 

0.6HNN  -  0.4XTHDCPD 

0.69 

0.29 

1.01 

1.38 

0.4HNN  -  0.6XTHDCPD 

0.75 

0.87 

1.02 

1.34 

0.2HNN  -  0. 8XTIIDCPD 

0.85 

0.66 

1.01 

1.20 

I 


Table  VI.  Test  ot  the  Rntti  -  Chau dh r i  model  for  the  NBC-DMO  bindry 
system 


i  ----- 

Mole  Fraction 

-W/R  from  nV 

-W/R  from 

i 

0.25HBC  -  0.75DM0 

11  +  7 

11  t  7 

50°  C 

1 

0.50NBC  -  0.50 DM0 

11  i  7 

11  +  7 

1 

0.75NBC  -  0 .25 DM0 

18  +  7 

24  1  7 

1 

1 

0.25NBC  -  0 .75 DM0 

21  t  7 

21  t  7 

o°c 

0.50NBC  -  0.50 DM0 

16  +  7 

17  +  7 

1 

0.75NBC  -  0.25DM0 

17  ±  7 

20  +  7 

1 

0.25NBC  -  0.75DM0 

19  +  ? 

20  +  7 

|  -50° C 

0.50NBC  -  0.50DM0 

20  +  7 

21  +  7 

1 

0.75 NBC.  0.25DM0 

24  t  7 

25  +  7 

I 

I 

r 

i 

( 

I 

1 
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Tableau,  Test  of  the  Katti  -  Chandhri  model  for  the 
system 


T°C 

Mole  Fraction 

-W/R  from  nV 

0.25CD  -  0.75 DM0 

128  +  6 

50°C 

0.50CD  -  0.50DM0 

140  +  6 

0.75CD  -  0.25DMO 

168  +  6 

0.25CD  -  0.75DM0 

159+6 

O 

o 

o 

0.50CD  -  0.50 DM0 

180  t  6 

0.75CD  -  0.25DM0 

209  t  6 

0.25CD  -  0.75DM0 

197  +  6 

•50°C 

0.50C0  -  0.50DM0 

226  +  6 

0.75CU  -  0.25DM0 

261  t  6 

CD-DMO  binary 

-W/R  from  n 
133  +  6 

145  +  6 

173 

164  +  6 

186  +  6 

215  t  6 
202  t  6 
232  +  6 

267  t  6 


»-t~ 
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Table  VIII.  Test  of  the  Katti  -  Chaudhri  model  for  the 
binary  system 


o 

T  c 


Mole  Fraction  -W/R  from  nV 


50°  C 


0GC 


-5Q°C 


0.201  INN 

-  0.80XTHDCPD 

22 

0.4  01  INN 

-  0. 60XTIIDCPD 

28 

0.60HNN 

-  0.40XTHDCPD 

18 

0.  HOI  INN 

-  0.20XTHDCPD 

12 

0.20HNN 

-  0.80XTHDCPD 

129 

0.40HNN 

-  0.60XTH0CPD 

143 

0.60HNN 

-  0.40XTHDCP0 

152 

0.80HNN 

-  0.20XTHDCPI1 

159 

0.20HNN 

-  0.80XTHDCPD 

476 

0.4  01  INN 

-  0. 60XTHDCPD 

512 

0.601  INN 

-  0.40XTHDCPD 

566 

0.801  INN 

-  0.20XTHDCPD 

643 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 


I1NN-XTHDCPD 

-W/R  from  n 
23  t  5 

31  i  5 

20  i  5 

16  t  5 

130  t  5 

146  +  5 

155  t  5 

162  t  5 

478  t  5 

514  i  5 

568  +  5 

645  1  5 
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Table^JXa.  VTF  Eq.  parameters  for  shear  viscosity  of  ND,  NBC,  CD,  NBB, 
DM0  and  their  mixtures  for  constant  B  parameter. 

In  n(P)  =  A  +  B/[T(K)  -  To] 


Molar  Composition 

T  range (°C) 

A 

B 

Jo. 

Std. 

Dev.  In 

ND 

-31 

to 

50 

-7.6864 

531.26 

98.9 

0.002 

NBC 

-57 

to 

50 

-7.5057 

II 

109.0 

0.006 

CD 

-55 

to 

50 

-6.9070 

II 

127.2 

0.033 

DM0 

-56 

to 

50 

-7.7613 

II 

97.5 

0.007 

NBB 

-55 

to 

50 

-7.6269 

II 

102.7 

0.014 

0.5ND  -  0.5NBC 

-40 

to 

50 

-7.5905 

II 

102.5 

0.002 

0.5ND  -  0.5CD 

-40 

to 

50 

-7.3976 

II 

108.1 

0.007 

0.5ND  -  0.5DM0 

-50 

to 

50 

-7.7245 

II 

98.3 

0.005 

0.5ND  -  0.5NBB 

-45 

to 

50 

-7.6391 

II 

96.1 

0.003 

0.5NBC  -  0.5CD 

-55 

to 

50 

-7.2497 

II 

116.6 

0.012 

0.75NBC  -  0.25DM0 

-54 

to 

50 

-7.5786 

II 

106.0 

0.007 

0.5NBC  -  0.5DM0 

-56 

to 

50 

-7.6397 

II 

103.1 

0.008 

0.25NBC  -  0.75DM0 

-55 

to 

50 

-7.7048 

II 

100.3 

0. 008 

0.5NBC  -  0.5NBB 

-53 

to 

50 

-7.5507 

II 

101.7 

0.005 

0.75CD  -  0.25DM0 

-54 

to 

50 

-7.2126 

II 

118.8 

0.014 

0.5CD  -  0.5DM0 

-55 

to 

50 

-7.4307 

II 

110.8 

0.006 

0.25CD  -  0.75DM0 

-53 

to 

50 

-7.6139 

II 

104.0 

0.004 

0.5CD  -  0.5NBB 

-54 

to 

50 

-7.3304 

II 

109.2 

0.006 

0.5NBB  -  0.5DM0 

-55 

to 

50 

-7.6821 

II 

96.1 

0.008 

1/3ND-1/3NBC-1/3CD 

-50 

to 

50 

-7.4243 

II 

109.1 

0.005 

7 


1/3ND-1 /3NBC-1/3NBB  -55 

to 

50 

-7.5935 

581.26 

99.4 

0.004 

1/3ND-1 /3CD-1/3NBB  -51 

to 

50 

-7.4667 

11 

103.6 

0.002 

1/3NBC-1/3CD-1/3NBB  -55 

to 

50 

-7.3916 

II 

108.4 

0.004 

0.25ND  -  0.25NBC  ,, 

-  0.25CL)  -  0.25NBB 

to 

50 

-7.4712 

II 

104.5 

0.003 

0.2MD  -  0.2NBC  -  0.2  „ 
CD  -  0.2NBB  -  0.2DM0  ~bb 

to 

50 

-7.5415 

II 

103.5 

0.005 

Hi) 


Table  IXb.  VTF  £q.  parameters  for  shear  viscosity  of  NBC-DMO  binary 
system  for  constant  B  parameter 


In  n  (P)  =  A  +  B/[T(K)  -  T.] 


Molar  Composition  T  range (°C)  A 


B  To 


Std . 
Dev.  In 


NBC 


-57  to  50  -7.3533  536.26  114.7  0.005 


0.75NBC  -  0.25DM0  -54  to  50  -7.4328 


111.9  0.005 


0.5NBC  -  0.50M0  -56  to  50  -7.4958 


109.1  0.005 


0.25NBC  -  0.75DM0  -55  to  50  -7.5639 


106.4  0.004 


DM0 


-56  to  50  -7.6226 


103.7  0.005 
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Table  IXc.  VTF  Cq.  parameters  for  shear  viscosity  of  CD-DMO 
binary  system  for  constant  B  parameter. 


In  n  (P)  =  A  +  B/[T(K)  -  To] 


Molar  Composition  T  range(0C)  A  B 


CD  -55 

0.75CD  -  0.25DM0  -54 

0.5CD  -  0.5DM0  -55 

0.25CD  -  0.75DM0  -53 

DM0  -56 


to  50  -7.1142  637.48 

to  50  -7.4041 

to  50  -7.6120 

to  50  -7.7881 

to  50  -7.9271 


Jo 

121.4 

112.4 

104.1 

96.9 

90.0 


Std. 
D_ev_1_  l  _n 

0.023 

0.006 

0.006 

0.008 

0.011 
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Tahiti  X.  Results  of  fits  to  master  viscosity  Eqs.  (43)  arid  (44). 

Units  of  Aj  ,  jj  Rare  so  chosen  to  siiear  viscosity  in  units 
of  poise. 


Mi xture 

com¬ 

ponent 

T  .  (K) 

A-j  values 
from  fit  to 
viscosity 

Eq.  (43) 

Std.  Dev. 
of  In  i, 
f  rom 

Eq-  (43) 

ai  values 
from  fit 
to  fluidity 
Eq.  (44) 

average  frac¬ 
tional  Std. 
Dev.  from 

Eq.  (44) 

ND 

96.63 

-7.715 

2255.02 

u 

NBC 

107.83 

-7.506 

1841.37 

Mixtures 

CD 

124.50 

-6.990 

0.066 

982.36 

0.050 

NBB 

97.60 

-7.638 

2086.22 

DM  0 

97.53 

-7.786 

2412.83 

B  = 

581.75 

B  = 

586. 

NBC 

114.65 

-7.359 

1562.81 

NBC -DM0 

DM0 

103.65 

-7.626 

0.007 

2043.95 

0.005 

B  = 

538.91 

B  = 

536. 

CD 

120.63 

-7.184 

1218.74 

CD-DMO 

DM0 

89.26 

-7.985 

0.064 

2807.66 

0.029 

B  = 

639.57 

B  = 

637. 

HNN 

152.27 

-7.081 

1224.54 

HNN- 

XTHDCPD 

XTH 

DCPD 

92.06 

-7.572 

0.056 

1932.18 

0.063 

B  =  830.34 


834.5 
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Table  XI .  Comparison  between  the  experimental  viscosity,  npxp,  and 
viscosity,  ncalc*  obtained  from  Eq.  (43)  using  tneM 
parameters  of  Table  X. 


100  (In  nexp  "  ncaic) 


CU)  Mixtures 

bu  C 

0°C 

-50  C 

NO 

5.0 

7.1 

9.6 

NBC 

2.0 

1.9 

5.1 

CD 

7.1 

18.3 

23.3 

DM0 

2.9 

1.5 

2.2 

HBB 

11.1 

10.2 

21.2 

0.5ND  -  0.5NBC 

2.1 

2.3 

2.6 

0.5ND  -  0.5CD 

-  6.4 

-  5.3 

-10.8 

0.5ND  -  0.5DM0 

4.2 

4.2 

7.2 

0.5ND  -  0.5NBB 

2.6 

1.3 

0.2 

0.5NBC  -  0.5CD 

-  1.5 

1.6 

1.4 

0.75NBC  -  0.25DM0 

-  2.2 

-  2.1 

-  3.1 

0.5NBC  -  0.5DM0 

1.7 

0.4 

2.2 

0.25NBC  -  0.75DM0 

-  1.4 

-  2.3 

-  2.6 

0.5NBC  -  0.5NBB 

1.3 

-  0.5 

-  2.0 

0.75CD  -  0.25DM0 

-  4.4 

-  3.6 

-  6.4 

0.5CD  -  0.5DM0 

-  5.1 

-  4.6 

-  5.7 

0.25C0  -  0.75DM0 

-  0.3 

-  2.8 

-  2.4 

0.5CD  -  0.5NBB 

-  3.4 

-  6.4 

-10.0 

0.5DM0  -  0.5NBB 

1.1 

-  0.9 

-  2.4 

8 


1/3  ND-1/3  NBC-1/3  CO 

-  3.7 

-  3.2 

-  5.7 

1/3  NO-1/3  NBC-1/3  NOB 

1.2 

-  0.3 

-  2.4 

1/3  NO-1/3  CD-1/3  NBB 

-  5.4 

-  7.5 

-13.2 

1/3  NBC-1/3  CO-1/3  NBB 

-  3.7 

-  5.0 

-  8.8 

0.25ND  -  0.25NBC  -  0.25CD 
-  0.25NBB 

-  3.7 

-  5.6 

-10.2 

0. 2ND  -  0.2NBC  -  0.2CD 
-  0.2DM0  -  0.2NBB 

-  2.8 

-  4.8 

-  7.1 

NBC-DMO 

5Q.°X 

Q°C 

-50°C 

NBC 

-  1.0 

-  0.9 

-  1.8 

0.75NBC  -  0.25DM0 

-  2.2 

-  2.1 

-  3.1 

U.5NBC  -  0.5DM0 

-  1.6 

-  2.3 

-  3.0 

0.25NBC  -  0.75DM0 

-  1.4 

-  2.3 

-  2.6 

DMO 

-  0.8 

-  1.2 

-  1.7 

CD- DMO 

50°C 

o°c 

-50C 

CD 

3.6 

9.7 

8.6 

0.75CD  -  0.25DM0 

-  4.4 

-  3.6 

-  6.4 

0.5C0  -  0.5DM0 

-  4.4 

-  6.4 

-  8.3 

0.25CD  -  0.75DM0 

-  0.3 

-  2.8 

-  2.4 

DMO 

7.2 

4.9 

7.5 
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I 

I 


HNN-XTJiDCPp 

| 

50°C 

o;'c 

-50°C 

! 

HNN 

• 

-  3.0 

-  4.6 

-13.3 

1  0.8HNN  -  0.2XTIIDCPD 

2.4 

2.3 

4.1 

■  0.6HNN  -  0.4XTI1DCPD 

4.0 

4.3 

9.8 

*  0.4 HNN  -  0.6XTHDCPD 

2.2 

3.3 

7.6 

|  0.2HNN  -  0.8XTHDCPD 

-  0.1 

0.0 

2.5 

XTHDCPD 

-  5.3 

-  5.2 

-  7.3 

t 

i 

! 

! 

I 

I 

l 

l 

{ 

f 

r 
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TabJe  XJI.  Comparison  between  the  experimental  viscosity,  neXp,  and 
viscosity,  nca^c,  calculated  via  Eq.  (44)  using 
the  parameters  of  Table  X. 


100  (In  neXp  -  In  ncaic) 


Cl0  Mixtures 

50°  C 

0°C 

-50°C 

ND 

3.8 

5.3 

6.8 

NBC 

1.3 

0.6 

2.6 

CD 

-  5.0 

5.4 

9.0 

DM0 

1.3 

-  0.6 

-  1.0 

NBB 

9.8 

8.3 

18.3 

0.5ND  -  0.5NBC 

1.6 

1.3 

0.5 

0.5ND  -  0.5CD 

-  4.8 

-  4.2 

-10.9 

0.5ND  -  0.5DM0 

2.8 

2.3 

4.3 

0.5ND  -  0.5NBB 

1.4 

-  0.5 

-  2.5 

0.5NBC  -  0.5CD 

-  3.1 

-  0.8 

-  2.1 

0.75NBC  -  0.25DM0 

-  0.5 

0.0 

-  0.2 

0.5NBC  -  0.5DM0 

1.5 

-  0.4 

0.3 

0.25NBC  -  0.75DM0 

0.2 

-  0.3 

0.2 

0.5NBC  -  0.5NBB 

0.4 

-  1.9 

-  4.4 

0.75CD  -  0.25DM0 

-  2.8 

-  1.6 

-  3.7 

0.5CD  -  0.5DM0 

-  2.2 

-  2.4 

-  4.6 

0.25CD  -  0.75DM0 

1.1 

-  1.0 

-  0.1 

0.5CD  -  0.5NBB 

-  3.2 

-  6.8 

-11.6 

0.5DM0  -  0.5NBB 

0.9 

-  2.7 

-  5.1 
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1/3  ND- 1/3  NBC-1/3  CO 

-  2.6 

-  2.7 

-  6.4 

1/3  ND- 1/3  NBC-1/3  NBB 

0.4 

-  1.7 

-  4.8 

1/3  ND-1/3  CD- 1/3  NBB 

-  4.0 

-  6.7 

-13.5 

1/3  NBC-1/3  CD-1/3  NBB 

-  3.4 

-  5.4 

-10.3 

0.25ND  -  0.25NBC  -  0.25CD 
-  0.25NBB 

-  2.8 

-  5.3 

-11.0 

0.2ND  -  0.2NBC  -  0.2CD 
-  0.2DM0  -  0.2NBB 

-  1.7 

-  4.2 

-  7.6 

NBC-DMO 

SOJLC 

0°C 

-50°C 

NBC 

-  0.1 

0.5 

0.4 

0.75HBC  -  0.25DM0 

-  0.5 

0.0 

-  0.2 

0.5NBC  -  0.5DM0 

0.3 

0.0 

0.0 

0.25NBC  -  0.75DM0 

0.2 

-  0.3 

0.2 

DMO 

0.2 

0.2 

0.4 

CD -DMO 

50°C 

-5JX°C 

CD 

-  3.0 

3.6 

3.2 

0.75CD  -  0.25DM0 

-  2.8 

-  1.6 

-  3.7 

0.5CD  -  0.5DM0 

-  0.9 

-  2.6 

-  3.8 

0.25CD  -  0.75DM0 

1.1 

-  1.0 

-  0.1 

DMO 

3.8 

1.8 

4.9 
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I 

( 


HNN-XTtIDCPD 

5o;x 

o'Jc 

(> 

-50  C 

HNN 

-  2.5 

-  5.1 

-16.1 

0.8HNN  -  0. 2XTHDCPD 

4.2 

3.3 

3.2 

0.6HNN  -  0.4XTHDCPD 

5.9 

5.6 

9.6 

0.4HNN  -  0.6XTHDCPD 

3.4 

3.9 

7.0 

0.2HNN  -  0.8XTHDCPD 

-22.5 

-  0.8 

-78.2 

XTHDCPD 

-  7.7 

-  8.1 

-11.0 

I 

I 

I 

! 

I 

! 

! 

I 

( 

l 

( 
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Table  XIII.  Comparison  of  molecular  diameters  of  C  hydrocarbons  calcu¬ 
lated  from  viscosity  with  those  estimated  from  density  and 
molecular  models 


Mol ecule 

mJa°1 

(M/^oNal^Aj 

L,  W,  H  ( A°) 
from  models 

a,>i  /  L 

ND 

9.5 

6.8 

12.8,  5.9,  5.6 

0.75 

NBC 

8.5 

6.6 

11.9,  6.3,  5.1 

0.71 

CD 

6.5 

6.3 

9.8,  7.5,  5.9 

0.67 

0M0 

9.9 

6.8 

13.4,  5.9,  5.1 

0.74 

NBB 

9.0 

6.3 

12.1  ,  7.0,  5.3 

0.74 

HNN 

7.1 

6.5 

9.3,  6.8,  6.8 

0.76 

XTHDCPD 

9.0 

6.2 

8.4,  6.9,  6.3 

1  .07 

nlcPi 


TI°CI 
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50 


CD -DM0 


I 

( 


I 

1  » 
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